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Although  tna  earliest  students  cf  hydrodynamics  v.ere  stronrly  interested 
in  such  natural  phenomena  as  fountains  and  water  waves,  propress  in  attaining 
a mathematical  description  of  these  phenomena  has  boon  relatively  slow* 

It  vas  rot,  in  fact,  until  the  nineteenth  century  that  the  mathematical 
framework  of  hydrodynamics  was  established  on  a firm  basis.  S-ince  then, 
knowledge  in  this  field  has  accumulated  — along  with  that  in  the  other 
major  sciences  — at  an  accelerating  rate.  Yet  there  remain  — as  in  the 
other  sciences  — important  questions  that  are  still  unanswered.  For  example, 
we  may  ask  what  are  the  precise  shape  and  velocity  of  a simple  surface 
wave  on  the  ocean?  Although  these  entities  can  be  approximated  to  a high 
degree  of  accuracy,  there  exist  as  yet  no  formulas  which  describe  them  exactly. 

The  basic  mathematical  difficulty  lies  in  the  treatment  of  the  free 
surface  boundary  of  the  fluid.  In  the  usual  problem  involving  a boundary, 
the  boundary  is  not  free,  but  is  instead  fixed  and  its  location  is  known. 

The  customary  approach  to  this  problem  then  involves  the  setting  up  of 
a mathematical  expression  which  states  that  the  fluid  mu3t  not  pass  through 
the  boundary.  However,  in  a hydrodynamical  problem  involving  a free  surface 
boundary,  the  character  of  the  boundary  value  problem  assumes  a special 
nature.  On  the  boundary  — the  free  surface  — one  prescribes  two  conditions, 
namely,  that  the  free  surface  is  not  violated  by  the  fluid,  and,  further, 
that  the  pressure  is  constant  there.  However,  until  we  have  solved  the 
problem,  we  do  not  know  the  shape  of  the  free  surface.  Thus  we  must  solve 
a boundary  value  problem  with  an  unknown  boundary. 

i 

Historically,  the  early  investigations  of  surface  waves  used  linearization 

i 

— a technique  in  which,  essentially,  the  free  surface  is  assumed  to  be 


flat  for  the  purpose  of  proscribing  boundary  conditions.  As  long  as  the 
waves  arc  very  small,  this  approximation  is  justified.  The  next  approach 
was  to  use  the  exact  equations,  but  to  neglect  all  terms  of  order  higher 
than  the  second.  The  applicability  was  thus  extended  to  include  larger  waves. 
Outstanding  among  the  contributions  to  the  theory  of  waves  of  finite  amplitude 
are  those  of  Stokes  [10]  and  Rayleigh  [8].  It  is  hardly  necessary  to  go 
into  the  history  of  thi3  theory,  since  excellent  treatments  can  be  found 
in  readily  available  works  such  as  Lamb  [3].  A less  technical  treatment 
is  given  by  Sverdrup,  et  al  [ll].  Rossby  [9]  has  shown  that  the  system  of 
deep  water  waves  of  finite  amplitude  given  by  Stokes  is  only  one  of  an  Infinite 
number  of  possible  irrotational  wave  systems.  Gemstner  f2]  dropped  the 
requirement  that  the  motion  be  irrotational  and  was  able  to  obtain  an  exact 
solution  to  the  remaining  equations  — those  of  continuity  and  constant 
surface  pressure  — which  showed  surface  waves  over  deep  water  to  have  a 
trochoidal  profile.  His  work  will  be  mentioned  again  in  Chapter  II. 

In  the  exact  theory  for  irrotational  waves  — the  theory  with  which 
this  paper  is  concerned  — Levi-Civita  [£]  has  shown  that  finite  waves  of 
permanent  form  do  exist.  It  might  appear,  especially  to  one  who  has 
witnessed  the  evidence,  that  an  existence  proof  for  ocean  waves  is  unnecessary. 
However,  this  is  not  the  case,  for  in  reality  it  is  only  in  the  mathematical 
treatment  that  true  permanent  waves  actually  exist.  We  neglect  viscosity 
and  hence  the  frictional  forces  that  in  nature  tend  to  cause  a wave,  in  time, 
to  lose  its  force.  Thus  in  deviating  slightly  from  physical  reality,  we 
make  it  necessary  to  prove  the  existence  of  a solution  to  our  problem.  The 
existence  proof  of  Levi-Civita  yielded  a fair  approximation  to  the  shape 
of  the  wave.  There  has  been  some  interest  in  the  solitary  wave  problem  also; 
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the  question  is  whether  or  not  a single  wave  of  permanent  form  can  exist 

on  the  ocean  inasmuch  as  a sudden  disturbance  usually  sets  up  a train  of 

waves  instead  of  a single  wave.  The  question  has  been  answered  in  the 

affirmative  by  Lavrentieff  [bj>  but  he  die  not  give  the  exact  shape  of  the 
# / 

wave»  This  problem  is  mentioned  because  of  the  fact  that  the  methods  of 
this  paper  could  be  applicable  to  it. 

A fresh  approach  to  free-surface  problems  in  hydrodynamics  has  been 
made  possible  recently  by  the  introduction  of  a new  mathematical  tool  [6]. 

This  new  procedure  allows  us  to  generate  free-surfaco  flows  at  will,  and 
it  forms  the  fundamental  concept  upon  which  this  paper  is  based.  The  method 
is  an  extension  to  the  study  of  non-gravity  flows  of  classical  investigations 
based  on  the  theory  of  functions  of  a complex  variable.  As  is  pointed  out 
in  most  textbooks  on  the  subject,  any  analytic  function  car.  be  interpreted 
as  the  solution  to  a two-dimensional  hydrodynamics!  flow  problem.  Our 
analysis  is  somewhat  reversed  from  normal  procedure,  however,  in  that  one 
must  first  decide  on  a function,  and  then  determine  what  hydrodynamics!  flow 
problem  it  solves.  Of  course,  more  than  pure  guesswork  is  involved.  Having 
a particular  flow  problem  in  rind,  one  proceeds,  ac  a r.  analytical  detective, 
to  uncover  certain  clues  which  will  enable  him  to  choose  the  proper  function. 
In  this  paper,  we  thus  construct  several  examples  f fr re-surface  flows  using 
a formula  suggested  by  Hans  Lewy  [6].  He  found  that  whereas  any  analytic 
function  represents  a hydrodynamical  flow,  if  v;e  restrict  the  functional  rela- 
tionship in  a certain  way,  the  flow  then  becomes  a free-surface  flow  in  a 
gravity  field. 

As  is  ofuen  the  case  with  mathematical  investigations  of  physical 
*/  This  theorem  has  also  been  proved  recent- ly  by  Friedrichs  and  Hyers. 
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phenomena,  one  must  make  some  plausible  assumptions  in  order  to  simplify 
the  mathematics.  In  this  paper  we  make  the  usual  assumptions  concerning 
the  fluid,  namely,  that  it  is  incompressible  and  is  without  viscosity. 
Further,  our  investigations  will  deal  with  phenomena  which  are  representable 
in  two  dimensions,  and  which  can  bo  generated  from  rest  — that  is,  which 
are  irrotational.  Bernoulli's  lav  ~iil  be  used  in  its  exact  form. 

Now  that  our  medium  is  defined  we  may  state  Lewy's  discovery  explicitly. 
He  found  that  if  Ac$)  is  an  arbitrary  analytic  function  of  ^ which  is  real 
on  the  real  ^-axis,  then  the  mapping 


(1) 


fd  A ( 3)*.2 

3? ; 


represents  a steady  free-surface  gravity  flow  in  the  z-plane,  the  free  surface 
being  the  image  of  a segment  of  the  real  ^-axis  along  which  the  inequality 


1 


(i/L)J 
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> 0 


holds.  In  Chapter  I,  this  relation  will  be  developed  in  a manner  different 
from  that  used  by  Lewy.  The  existence  of  /)  for  such  a free-surface  flow 
and  the  proof  <~f  the  analyticity  of  the  free  surface  will  also  be  given. 

It  should  be  pointed  out  that  this  work  was  done  by  Lewy  in  a more  general 
manner  than  that  presented  hem.  It  is  hoped  that  the  loss  in  generality 
is  compensated  for  by  a gain  in  simplicity. 

The  major  contribution  of  this  paper  lies  in  the  application  of  Lewy's 
theory.  Up  until  +his  time  there  have  been  no  important  free-surface  flows 
in  a gravity  field  represented  in  closed  form.  The  solutions  have,  in  the 
past,  been  either  approximations  (linearized  theory)  or  the  first  few  terms 
of  an  infinite  series  solution.  With  Lewy's  method  one  obtains  a solution 


■which  is  in  closed  form*  The  major  difficulty  with  the  method  lies  in  trying 
to  pick  the  few  significant  functions  /(  ( ) out  of  the  infinite  number 
of  possible  functions*  The  functions  wuich  will  yield  the  solutions  to 
three  outstanding  problems  of  surface  waves  remain  as  yet  unknown.  These 
are  the  problems  of  periodic  waves  over  a flat  bottom  of  finite  depth, 
periodic  waves  over  an  ocean  of  infinite  depth,  and  the  solitary  wave. 

It  is  hoped  that  this  paper  will  be  of  some  aid  in  the  quest  for  the 
solutions  to  these  problems • 

In  Chapter  II  an  example  of  a free-surface  wave  will  be  given.  It  is 
of  trochoidal  profile  and  travels  over  a smooth,  but  not  a flat,  bottom. 

Chapter  III  contains  a limiting  case  of  the  flow  of  Chapter  II  where  the  trochoid 
becomes  a cycloid.  It  is  shown  that  when  this  occurs,  the  character  of  the 
flow  changes  completely  in  that  the  periodic  nature  of  the  flow  is  destroyed. 

Some  general  properties  of  deep-water  waves  are  studied  in  Chapter  IV  with 
respect  tc  their  affect  on  the  choice  of  the  arbitrary  function/?  (3). 

An  example  of  surface  waves  growing  out  of  these  considerations  is  given. 

In  Chapter  V we  leave  the  realm  of  surface  waves  and  examine  two  flows 
of  a non-periodic  nature.  One  is  a fountain-like  flow  which  yields  many 
interpretations.  It  is  especially  interesting  because  it  provides  a first 
example  of  how  a stagnation  point  on  the  surface  need  not  necessarily  occur 
at  an  angle  of  2 TT/3 » The  second  example  might  be  interpreted  as  the  skimming 
action  of  a vertical  board  being  drawn  rapidly  along  the  surface  of  a deep 
body  of  fluid. 

It  should  be  mentioned  that  the  attack  on  the  problems  of  free-surface 
phenomena  received  impetus  during  the  war  because  of  its  military  importance, 
the  major  interest  centering  around  ocean  waves.  In  order  to  prepare 
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adequately  for  an  amphibious  assault,  one  needed  a means  of  forecasting  wind 
waves  and  swell.  A part,  of  this  study  required  a knowledge  of  shapes  and 
velocities  of  the  waves  when,  having  betr  created  by  storm,  they  moved  on 
through  the  undisturbed  sea.  A knowledge  of  the  shape  and  velocity  of  ocean 
waves  is  also  important  to  those  making  a study  of  ship  stabilization, 
since  the  forces  on  the  hull  of  a ship  at  sea  are  dependent  on  them.  Under 
certain  conditions,  one  can  see  ocean  waves  quite  clearly  from  the  air. 

To  the  airman,  the  exact  shape  of  the  waves  is  not  too  important,  but  the 
velocities  might  well  concern  him.  A ship  at  sea  forms  directly  behind 
it  a train  of  waves  called  the  transverse  wake  and  these  waves  move  with 
the  same  velocity  as  the  ship.  From  the  air  one  can  measure  the  wavelength 
of  this  transverse  wake  and  consequently,  knowing  the  velocity  of  a wave 
as  a function  of  its  wavelength,  one  is  able  to  determine  the  speed  of 
the  ship.  Also,  the  knowledge  of  wavelength  as  a function  of  velocity  and 
depth  can  be  applied  to  determine,  from  aerial  photographs,  the  depth  of 
a bay,  the  slope  of  a beach,  and  other  geographical  data. 

I would  like  to  express  my  great  appreciation  for  the  help  given  by 
Professor  Hans  Lewy  in  my  early  work  on  the  applications  of  his  theory. 

And,  of  course,  the  guidance  and  inspiration  received  from  ray  advd'or 
Professor  Paul  Garabedian  is  without  measure. 
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G'uapter  I 

As  wa s mentioned  in  the  introduction,  we  are  investigating  steady, 
two-dimensional,  irrotational , free-surface  flows  of  an  incompressible, 
inviscid  fluid  in  a gravity  field.  From  the  irrotationality,  we  know  that 
there  exists  a velocity  potential  p such  that  the  velocity,  q , is  the 
gradient  of  p . For  steady  incompressible  flows  the  equation  of 

continuity  states  that  the  divergence  of  "q^is  zero;  therefore  <p  satisfies 

. . . . 2 „ 

lapuaces's  equation  v 'P  = 0. 

Bernoulli's  la*. , which,  in  our  case,  takes  the  form 
(2) 


Finally,  the  dynamical  equations  yield 


2 P 

1/2  q ♦ gy  ♦ ^ = constant. 


f 

— > 

where  q , without  the  arrow,  is  the  magnitude  of  the  velocity  q , g is  the 
acceleration  due  to  gravity,  y is  the  vertical  coordinate  (positive  upwards), 

P is  the  pressure  at  any  point  of  the  fluid,  and  ^ is  the  density  of  the 
fluid.  If  we  introduce  the  stream  function  we  may  describe  a two- 
dimensional  flow  by  the  complex  potential 

^(z)  (z)  ♦ i^y/(z). 

Then 

ar*  ’ f x * 1 % • fx  - 1 f y * "i  - 1 % • 

where  q ^ and  are  the  x and  y components  of  velocity,  respectively,  whence 

— 

q 

the  bar  denoting  the  complex  conjugate.  The  above  development  is  familiar 
from  textbooks  on  hydrodynamics  or  complex  variables  and  consequently  we 
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omit  further  detail  nere  and  proceed  to  develop  Lewy's  theory. 

Fnysicaiiy  a stream  lire  is  a line  across  which  no  fluid  will  flow  — 
thus  any  stream  line  is  a possible  boundary  for  the  flow.  If  we  want  one 
of  our  stream  lines  to  be  a free  surface  boundary,  then  we  must,  in  addition, 
require  that  the  pressure  be  constant  all  along  this  stream  line.  Conversely, 
if  along  some  segment  of  a stream  line  the  pressure  turns  out  to  be  constant, 
then  it  is  a possible  free  surface.  For  the  best  physical  example  to 
illustrate  this  constant -pressure  criterion,  we  again  turn  to  the  ocean. 

The  constant  pressure  along  the  free  surface  is  ,just  equal  to  the  atmospheric 
pressure  of  the  air.  Since  the  inertia  of  air  is  very  small  compared  to 
that  of  water,  the  surface  of  the  water  may  easily  displace  the  air  and  so 
is  not  impeded  by  the  air  during  a change  in  shape.  That  is  to  say,  the 
surface  is  free  to  change  shape.  This  is  an  important  factor  in  progessive 
waves , for  if  they  were  required  to  do  much  work  in  displacing  the  air,  they 
would  soon  lose  their  force.  Of  course,  the  constant- pressure  criterion  is 
again  a mathematical  idealization,  for  air  does  have  its  effect  on  surface 
waves,  especially  when  the  air  is  itself  in  motion.  Indeed,  most  ocean  waves 
owe  their  existence  to  the  forces  created  by  air  in  rapid  motion. 

In  any  flow  represented  by  a complex  potential,  a stream  line  is 
characterized  as  the  image  of  a line  Im( 5 ) = constant.  One  of  the  stream 
lines  in  the  z-plane  is  hence  the  image  of  the  line  Tm(^)  * 0,  and  it  is 
this  line  that  we  will  select  for  the  free  surface.  Tie  now  turn  our  attention 
entirely  to  this  free-surface  stream  line.  It  may  be  that  only  part  of  the 
real  J’-axis  will  satisfy  the  free-surface  criterion  of  constant  pressure. 

At  any  rate,  whether  it  be  all  or  only  a part  of  the  ^-axis  that  is  the  image 
of  the  free  surface,  we  will  call  it  the  free-surface  segment  of  the  real^^-axis. 
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Along  this  free-surface  segment,  or,  more  properly,  along  its  image  in  the 
z-plar.e,  the  pressure  P is  constant*  Since  the  fluid  is  incompressible,  the 


UCiiOlbjr  j JD  d bUJIO  uauv  biawuf.  uv 


ghoui  the  fluid  and,  therefore,  in  particular,  on 


the  free  surface.  We  may,  there  fore,  write  Bernoulli's  t-’w  in  the  form 

2 P 

(3)  1/2  Q ♦ gy  = y ♦ constant, 

where  everything  on  the  right-hand  side  is  constant.  We  are  at  liberty  to 
chocse  the  location  of  the  origin  in  the  z-plane  and  can  do  so  in  such  a way 
that  the  constant  gk  obtained  from  replacing  y by  some  y*  ♦ k just  cancels 
the  two  constants  on  the  right-hand  side  of  equation  (3).  This  means,  in 
particular,  that  the  zero-velocity  level  of  the  fluid  coincides  with  the  real 
z-axis  when  we  are  on  a free  surface.  Notice  that  if  we  are  not  on  a free 
surface,  there  is  no  reason  to  expect  zero  velocity  on  tne  real  z-axis,  since 
then  the  pressure  P is  not  necessarily  constant.  A saving  in  effort  is 
achieved  if  we  further  set  the  acceleration  of  gravity  equ  1 to  1.  Physically 
this  means  we  are  choosing  our  unit  of  length  to  be  equal  to  approximately 
32  feet,  provided  time  is  still  measured  in  seconds.  Thus,  on  the  free  surface, 
Bernoulli's  law  takes  the  simple  form 


^ I 2 


♦ v = 0 


This  can  be  inverted  and  rewritten  as 


dz  dz  , 1 

^ 2y 


dj  2y 

This  equation  is  true  only  on  the  free  surface,  with  dz  and  d j measured  along 
the  free  surface.  But  the  free  surface  is  the  image  of  Im(^)  « 0,  and 
therefore  we  may  replace  d~^  by  d , since  both  are  equal  to  d and  are 
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real  there.  Hence  equation  (5)  becomes 


dz  dz 

d T dT 

S>  -J 


2y 


Since  z - z = 2iy,  v:e  have  z = z - 2iy,  so  that 


uz 


dz  


2i& 

<i£  ’ 


dy 


Q-  has  meaning,  since  the  free  surface  can  be  represented 

ap 


wnere 

parametrically  in  the  fcrra  z = z((p)  = x((p)  * iy((p  ) . We  note  that  from 
equation  (U),  the  inversion  z * z( j)  is  always  possible  if  the  free  surface 
is  bounded  away  from  the  real  z-axis.  If  we  let  a prime  1 designate  differentia- 
tion with  respect  to  we  have  finally 

1 


(6) 


z'  (z!  - 2iy! ) ♦ 


w 


o . 


We  now  come  to  the  essence  of  Lewy's  idea.  Equation  t6)  is  an  equation, 
involving  a real  quantity  y,  that  must  hold  on  the  free  surface.  But 
suppose  we  consider  it  as  an  equation  in  the  complex  domain  involving  a 
complex  function  y.  Then  its  solution  would  be  an  analytic  function 
(under  proper  restrictions)  and  hence  would  represent  a hydrodynamical  flow. 
If  this  complex  quantity  became  real  along  the  stream  line  Im(^)  *=  0,  then 
by  virtue  of  the  derivation  of  equation  (6),  that  stream  line  would  be  a 
free  stream  line.  To  accomplish  this,  we  replace  -y  Dy  cm  analytic  function 
of  5,  which  is  real  and  positive  along  some  segment  of  the  real 

■^-axiSj  If  we  insert  this  function  into  equation  (6)  we  have 

(?)  z'  (*»  ♦ 2i  J')  - « 0 . 

Since  A (^)  is  analytic,  equation  (7)  is  a non-linear  differential  equation 


! 


JJ 

which  connects  z with  ^ in  an  analytic  way.  Consequently,  it  determines  a 
hydrodynamical  flow.  Furthermore,  the  image  of  a segment  of  the  real^-axis 
is  a stream  line  along  which  pressure  is  constant  by  virtue  of  the  derivation 
of  equation  (7),  provided  /l  = -y  on  the  stream  line  — a fact  we  establish 
immediately.  Hence  this  stream  line  is  a free-surface  stream  line. 

The  differential  equation  (7)  may  be  solved  for  z in  terms  of  an  integral, 

(8)  z - -i/l  * f \~rr  - d ^ * 

\J  ' — 

It  is  clear  from  this  equation  that  A (^)  is  equal  to  -y  on  the  free 
surface,  provided  the  integral  itself  is  real  on  the  free  surface. 

Consequently^ we  must  require 

(u  -It-  - *'2-0 

along  the  free-surface  segment.  We  call  this  the  free-surface  condition. 

We  now  have  a method  of  constructing  free-surface  flows.  Using 
equation  (8)  we  may  generate  a free-surface  flow  by  inserting  any  function 
^ (^ ) which  is  1)  analytic,  2)  greater  than  zero  (hence  real)  on  some 
segment  of  the  real  t^-axis,  and  3)  such  that  1/2/^  > on  this  segment. 

Recall  that  in  the  classical  theory,  the  arbitrary  analytic  function  itself 
determines  the  flow,  whereas  here  the  arbitrary  analytic  function  appears 
as  a parameter  and  as  a consequence  exerts  a more  indirect  influence  on  the 
resulting  flow. 

It  is  also  possible  to  make  a converse  statement.  For  any  free-surface 
flow,  there  exists  a unique  function  ) which  will  generate  the  flow  in 

tne  above  sense.  Furthermore,  any  free-surface  flow  may  be  continued  across 
the  free  surface  and  as  a consequence,  the  free  surface  must  be  an  analytic 


wmmmr*?  ■ - 


- 12- 


curve.  Before  goinp  ahead  with  the  proof  of  these  statements,  we  will 
demonstrate  the  necessary  existence  and  uniqueness  theorems  for  first-order 
ordinary  differential  equations  in  the  complex  domain. 

The  existence  and  uniqueness  theorems  needed  here  differ  slightly  from 
those  usually  given  in  the  textbooks.  Often  the  initial  values  are  taken 
the  centers  of  the  domains  under  consideration.  v»e  will  need  to  assign 
an  initial  value  to  the  independent  variable  which  is  on  tne  boundary  of  the 
domain.  Further,  the  usual  proof  just  shows  that  the  solution  exists  and 
is  analytic  within  the  domain  under  consideration.  Vfe  will  need  both  existence 
and  uniqueness  on  the  boundary  of  the  domain,  and  we  will  prove  that  the 
solution  is  also  continuous  on  the  boundary. 

Consider  the  differential  equation 

(10)  * f(^  ) 

together  'with  the  following  domains.  Let  be  an  open  circular  region  of 
the  ^ -plane  of  radius  r and  center  /?Q.  Let  D2  be  an  open  simply-connected 
region  of  the  T^-plane  which  may  be  of  any  configuration,  provided  that  any 
two  points  belonging  to  can  be  joined  by  a path  of  finite  length  contained 
in  the  closure  of  D?.  Choose  the  radius  r of  and  the  domain  D?  so  that 
f(/)  is  regular  and  single  valued  for  both  variables  ranging  within 

their  respective  domains,  is  continuous  in  both  variables  within  and  on  the 
boundaries  of  these  domains,  and  satisfies  the  Lipshitz  condition 

f (_$,/?!)  ~ f(5,/?2)  <K^1-/?2 

for  A and  ^2  or  on  boundary  of  and  ^belonging  to  the  closure 

of  Dg.  Further,  let  M be  the  maximum  value  of  f(^>/|)  for  y and  A belonging 


|. 


(v  »rv>*  i -■« 


■ I 
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to  the  closures  c’"  their  respective  domains . The  existence  of  the  maximum 
is  assured  by  the  continuity  of  the  function  f(  gj).  Now  let  belong 

to  the  closure  of  the  domain  DQ,  -where  we  specifically  allow  values  on  the 
boundary  of  D2„  If  we  denote  the  greatest  length  of  path  necessary-'  to  join 
^ to  any  point  j that  belongs  to  the  closure  of  Dg  by  S,  and  if  S < I ' 
then  we  may  proceed  with  the  analysis.  If,  however.  S > ^,then  we  alter  D2 
by  removing  those  parts  farthest  away  (by  path)  from^Q  until  the  inequality 
S < £ holds.  Call  this  altered  domain  the  domain  D,  and  set  D * D,  when 
no  alterations  are  necessary.  Then  the  differential  equation  (10)  admits 


a unique  solution 


with  /' 


' >•  > 


which  is  analytic  within  D and 


continuous  in  the  closure  of  B. 


Suppose  for  a moment  that  a solution  to  the  differential  equation  is  known 
which  has  as  initial  value  for  \ = This  solution  would  satisfy  the 


relation 


A(^)  - ?0  ♦j.  f(t,  A (t))dt. 


If  the  function  is  unknown,  then  this  relation  becomes  an  integral 

equation  which  may  be  solved  by  the  method  of  successive  approximations  in 
the  following  manner. 

Consider  the  sequence  of  functions 


^Although  it  is  not  really  necessary,  we  will  assume  for  sake  of  clarity 
that  a geodesic  exists  for  all  pairs  of  points  belonging  to  the  closure  of 
D2,  We  say  geodesic  (path  of  shortest  length)  because  it  is  possible  that 

two  points  could  not  be  joined  by  a straight  line  lying  in  the  closure  of 
the  domain  B2*  In  that  c'^e  parts  of  the  boundary  must  be  used  together 

with  straignt-line  segments. 


I 


• c 

A,<5>  - * j”  f(t,^)dt, 

A„d>’A0*J  r(t,  A^Jdt, 


where  each  of  the  integrals  is  a line  integral  from  to  some  point  ^ in 

D along  a path  each  point  of  which  belongs  to  the  closure  of  D.  Since  /? 

is  a constant  within  D, , f(t,  ) is  a regular,  single- valued  function  of  t 

• 1 o 

for  t within  D}  and  is  continuous  in  the  closed  domain  D.  3y  Cauchy's  integral 
theorem  the  value  of  ^ is  independent  of  the  path  of  integration  and, 
further,  by  Morera's  theorem,  /\ ^ is  regular  and  single  valued  within  D.  The 
continuity  of  ^ for  all  ^belonging  to  the  closure  of  D follows  from  the 
fact  that  it  is  an  integral.  Note  that  the  path  of  integration  may  include  part 
of  the  boundary  of  D,  for  Cauchy's  theorem  requires  only  analyticity  within 
the  contour  and  continuity  on  the  contour.  The  endpoint  may  be  on  the  boundary, 
and  the  value  of  /?.,  on  the  boundary  is  quite  independent  of  our  path;  that 
is,  we  may  reach  this  boundary  point  by  integration  along  part  of  the  boundary 
or  by  approaching  it  from  within  D. 

Turning  our  attention  to  we  see  that  we  will  be  able  to  attribute 

the  same  properties  to  as  we  have  to  as  soon  as  we  have  shown  f(t, 
to  be  regular  and  single  valued  within  D and  continuous  for  ^ belonging  to 
the  closure  of  D.  To  do  this  we  must  first  show  that  is  within  D^.  Now 


V*. 


)dt. 


path  of  integration.  Consequently, 


-l£- 


where  the  value  on  the  left  is  independent  of  the 

for  purposes  of  obtaining  a good  estimate,  we  may  choose  the  geodesic  joining 

* * 

to  as  our  path  of  integration. 


\ - ^ o i ( 

' ^4 


-<u  J; 


•j 
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dt 


dt 


< MS 


j* 

But  by  construction,  S < r;  so  that 

— M 

ana  is  in  D^.  #e  know  that.  f(^.A)  is  regular,  continuous,  and  single 
valued  in  its  second  argument,  provided  the  second  argument  is  within  D^. 

But  ^ ia  within  D^;  consequently,  f( $ , ^ ) is  a regular  single-valued 
function  of  ^within  D and  is  continuous  in  the  closure  of  the  domain  D. 

From  this  we  may  conclude  that  /? ^ is  a regular  single-valued  function  of 
$ within  D and  a continuous  function  of  J^within  and  on  the  boundary  of  D. 

Suppose  now  that  /i  ^ is  within  and  is  a regular,  single-valued  funct 
of  ^within  D,  and  is  continuous  in  the  closure  of  D.  Then 


i 


y 

~b 


dt 


< MS  < r, 


so  that  /{  also  is  within  D^.  Further, 

K = ^o  *5  f(t*^n-l}  dt  ’ 


ion 


/ 


so  that  /\  is  also  regular  and  single  valued  within  D and  continuous  in  the 

closure  of  D,  Thus  the  induction  is  established. 

o 

We  now  show  that  these  /\  converge  uniformly  and  absolutely  to  a limit 
function  /?(<).  In  the  remainder  of  this  proof  we  will  always  agree  to  take 
the  geodesic  path  when  integrating  to  a point  ^ . It  is  clear  that 

h -.1  I ^ I „ II  I 


/\-  A < f f(t,/)0)  « , 


6 

<M\  ds  ® Ms  , 

where  s is  the  distance  from  to  1?  along  the  geodesic  path.  Notice  that 
3 is  a function  of  J and  is  not  to  be  confused  with  S,  which  is  a constant. 
If  D is  convex,  then  s - J0j  . Now 

A2  - < r kt,  dt , 

J -e 


p 

l 


and  upon  applying  the  Lipshitz  condition  we  have 

u -,1  pi  i ^ I i 


?i-K 

I V a 


But  we  have  just  shown  that  < Ms , so 


! ff*  * 
1' 

1 I 

j f 


< 


MK  \ s ds  * MK  -S-r  . 


Suppose  now  that 


^n-1-  'i-2|  <MKn"2  ? 


chut: 


i § 

ft 


1 
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We  have 
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»S 

n 

' n-1 

\ 

- " | ^n-l  - Ki 


dt 


1 


dt 


t — TTT~  \ sn"*1  ds  ” MK11"1  “T 

( n~l) ! J nj 

5. 


IF, 

' 

1TV  , 

$ i 
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ft 

K 

( 

ll 

1 1 
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We  may  now  replace  s by  S sc  that  we  have  a uniform  bound  for  the  absolute 

values  of  the  differences  /{  - fl  , . Consequently  the  3um 

n n— i 


\ A,<^>  ' ■ 


by  the  Weierstrass  M-test,  converges  absolutely  and  uniformly  in  the  closure 
of  the  domain  D , Therefore,  An—y2(l)  and  }( 5 ) i3  regular  and  single  valued 
within  D and  continuous  within  and  on  the  boundary  of  D,  because  a uniformly 
convergent  sequence  of  continuous  or  regular  functions  converge  to  a continuous 
or  regular  function. 

Consider  the  equation 


.5 


/?n(!^)  “ P-Q  * ^ f(t,  ^ n_l.(t)  )dt  . 


S 


If  we  take  the  limit  as  n — ±cn  of  both  sides  we  obtain 


s 


lim  -^n(3)  * /?0  ♦ iira  ^ f(t,  /^^(t^dt  , 

n— >oo  ^ - 


n — » oo 


X 


! 


/ 


This  last  step  is  legitimate  by  virtue  of  the  uniform  convergence  of  the 

functions  P . 

n 

If  «e  consider  the  Lipshitz  condition 


< K - /?n  <K£n, 


we  see  that  the  functions  f(3,  /10),  f ( J , /l^),  f(o  , /C,)  ...  form  a uniformly 
convergent  .sequence  of  regular  single -valued  functions  within  D,  continuous 
in  the  closure  of  O.  Consequently  the  limit  function,  is  also 

regular  and  single  valued  within  I*  and  continuous  in  the  closure  of  D.  Under 
these  conditions  we  may  differentiate  the  indefinite  integral  (ll)  to  obtain 

Thus  P\  ) satisfies  the  differential  equation  and  from  (11)  it  is  clear  that 

A < V - 

"Finally,  we  prove  the  uniqueness  of  /?  . be  a function  which 

assumes  the  value  P at  and  which  satisfies  equation  (10),  and  hence 
equation  (ll).  Suppose  that  is  regular,  single  valued,  and  continuous 
in  a connected  domain  D*  which  includes  those  regions  of  D in  which  the 


inequality 


holds.  At  least  a small  region  of  D near  is  included  by  reasons  of 


-19- 


continuity.  If  ^ is  now  restricted  to  l>w,  we  have 

, »S 

J\  - A - ( rf(t.y\)  - f(t.^)]dt  . 
0 S 


whence 


A - /l|<  r |f(t,A)  - f(t, 

vj  y? 

< K A | A-  A\  |dt|, 

o 


o5 

< Kr  ) ds  ■ Krs 
J o 


where  s has  the  sane  meaning  as  before.  But  again 

/1-A\  < | dt , 

•-Vi 


< K J A -/i||dt|. 


We  may  now  use  the  estimate  above  to  obtain 


>f  2 32 

s ds  - IT  r !y 


It  is  quite  clear  that  the  induction  can  be  established  so  that  we  obtain, 
finally, 
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n 0n 

- .,n  ^ s r.n  S 

— nl  — nl 


But  the  left-hand  side  is  independent  of  n;  consequently,  A.”  A in  D*. 

But  if  3-A  in  D»>then  necessarily  /\  - ./"Yin  D.  This  completes  the  proof. 

We  now  have  the  equipment  to  prove  that  a unique  function  A (^)  exists 
for  any  free-surface  flow,  and  that  the  free  surface  is  an  analytic  curve-r 
Our  basic  tool  is  the  differential  equation  (7).  We  consider  as  given  a complex 
potential~5  (z)  = (p  ♦ i"^which  represents  a free-surface  flow,  namely, 
which  is  analytic  in  some  region  of  the  z-plane  bounded  by  an  arc  lying  in 
the  lower  half-plane  whose  image  is  a segment  of  the  real  -axis.  Along 
this  segment  we  assume  that  the  equation 

1d5|2 


1/2 


d z 


y - o 


holds „ This  is  tantamount  to  saying  „hat  the  pressure  is  constant  along  the 
arc.  This  exceptional  arc  will  be  called,  and  is,  the  free  surface. 

For  the  proof,  it  suffices  to  consider  a small  arc  of  free  surface  which 
is  bounded  in  depth  and  which  is  bounded  away  from  the  real  z-axis.  That 
is  to  say,  on  this  free-surface  segment,  M > y > f.  > 0.  It  is  then  clear 

from  Bernoulli's  law  that  z'  = is  also  bounded  and  bounded  away  from 

■at 

zero.  We  will  further  assume  that  z'  is  continuous  on  the  free  surface. 

This  means  that  the  free  surface  must,  in  particular,  have  a continuously 
turning  tangent.  It  is  given  that  z'  is  analytic  within  some  semi-neighborhood 


—The  method  used  in  proving  this  theorem  has  been  presented  in  more  generality 
by  Lewy  in  A theory  of  terminals  and  the  reflection  laws  of  partial  differential 
equations,  Technical"Tte  port  No.  Y,  Contract  Nonr-22^ (il) (NR-O/4I-O86) , Office  of 
Naval  Research,  Applied  rat-hematics  and  Statistics  Laboratory,  Stanford  Uni- 
versity, California,  Aug.  8,  1V5'2. 
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of  the  free  surface. 

We  first  find  the  /{  function  that  goes  with  this  flow.  Equation  (7), 


(7) 


(z*  ♦ 2i  A1)  - - 0 , 


may  be  considered  as  a differential  equation  for  an  unknown  /?(^)  with  Known 
z'.  Solving  for  A ! we  have 

(12)  A'  * d^~  m Ti  'K  2z'^  “ z')* 

For  initial  values,  we  will  take  a point  on  the  free  surface.  Let  /?  =•  - y , 

where  y^  is  the  imaginary  part  of  a point  zq  on  the  free  surface  which  is  the 

a 

image  of  j The  region  is  some  lower  semi- neighborhood  of  the  image  of 

+ha  2Tc  undsr  cor.sidsrstion  z*  is  ^ continucuc  ^und 

single  valued.  It  is  clear  that  all  the  hypotheses  about  the  function 

f(-5 , 2)  - -|j;  (-  — — z * ) are  satisfied  if  we  take  for  any  circle 

which  excludes  the  origin  of  the  A -plane  and  for  some  lower  semi-neighborhood 

of  the  real  ^-axis  where,  by  reason  of  continuity,  z'  remains  bounded 

away  from  zero  and  infinity  as  it  is  on  ths  free-  surface  segment.  Thus 

n 

there  exists  a unique  function  which  satisfies  equation  (12)  and 

attains  the  value  - y at  . 

o -n 

To  see  that  this  A ("^)  becomes  -y  all  along  the  free-surface  segment 
included  in  D and  not  just  at  the  point  zq,  we  use  the  uniqueness  of  the 
solution.  It  suffices  to  show  that  -y  is  a solution  to  equation  (7)  along 
the  free-surface  segment.  Consider  Bernoulli's  law  in  the  form 

*1  ♦ -A_  - o . 

d£  ?y 

Since  we  are  only  interested  in  staying  on  the  free  surface,  = d ^ ; 
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« dz 

hence  - d ■{  , and  >ve  may  replace  — — by  -~j  . Further,  z-z  = 2iy,  and 

d £ a 5 

on  the  free  surface  we  have  assumed  that  the  real  and  imaginary  parts  of  z 
(now  functions  of  the  real  variable  (p  alone)  have  continuous  derivatives,  so 


that  z'  = z ’ - 2iy'  has  meaning,  where  prime  designates 


Bernoulli's  law  takes  the  form 


3J"  T 


. Thus 


z'  (z'  - 2iy«)  4 _±_  = 0 , 


which,  upon  replacing  y by  -y,  is  just  equation  (7).  This  verifies  that  -y 
is  a solution  of  equation  (7)  and  therefore  that  /{  (^)  and  -y  coincide  on 
the  free  surface.  Incidentally,  this  shows  that  is  independent  of  the 

starting  point  zq. 

As  a function  of  y /?  is  real  on  the  real  ^ -axis  and  hence  can  be 
continued  into  the  upper  half-plane  by  Schwarz's  reflection  principle.  This 
defines  /{  (^)  as  an  analytic  function  above  and,  in  particular,  on  the  real 
^-axis.  We  then  return  to  equation  (7)  and  consider  it  as  a differential 
equation  in  an  unknown  z(-^)  and  a known  /(  (^)  above  the  free  surface. 

As  such  it  is  solvable  in  the  form  of  equation  (8),  which  thus  serves  to  deter- 
mine z above  the  free  surface.  Because  this  new  z(j^)  defined  above  the 
free  surface  agrees  along  the  free  surface  with  the  original  z(-^)  of  the  flow, 
it  is,  b;  definition,  the  analytic  continuation  of  the  original  z(^).  But 
being  an  analytic  function,  determines  a hydrod.’Tiamical  flow. 
Consequently^ we  have  extended  the  flow  across  the  free  surface  boundary. 

The  continued  function  del)  is  analytic  on  the  free-surface  segment,  so 
byMorera's  theorem  and  equation  (8),  z(^)  is  also  analytic  on  the  free 
surface.  Since  the  free  surface  is  the  image  of  Im(^)  - 0,  it  is  represented 


in  an  analytic  way  by 


z = z(  (p)  * x(<p)  ♦ iy(p)  , 


I 


~~ 7 


and  is  hence  an  analytic  curve. 

This  proof  breaks  down  completely  if  we  allow  the  free  surface  to  reach 
the  line  y * 0,  that  is,  if  we  allow  a stagnation  point  to  occur  on  the 


surface,  for  then 


2=0 


and  the  Lipshitz  condition  cannot  be  satisfied.  As 


Stokes  pointed  out,  the  free  surface  may  then  have  a cusp  of  angle  — ^ — at 
this  point.  This  can  be  seen  as  follows:  Let  the  point  in  question  be  the 
origin  in  the  f-.  R-.  and  z-planes.  If  the  slope  of  the  free  surface  near 
the  stagnation  point,  is  bounded  away  from  zero  and  infinity,  then  x and  y 
will  be  of  the  same  order  there.  Consequently,  z will  be  of  the  same  order 
as  that  is, 

z and  2”^  h2(^>)  > 


where  R^  and  R,,  are  regular  and  bounded  away  from  zero  in  the  partial 
neighborhood  of  the  origin  within  the  region  of  flow.  For  the  differential 
equation  (7)  to  be  satisfied,  the  order  of  the  first  term  must  be  the  same 
as  the  order  of  the  last.  If  we  substitute  the  values  of  z and  2 given 
above,  we  obtain 

where  R^,  Ri  5 and  have  the  same  properties  as  R,  and  R».  The  order  of 
i 4 b l d 

the  first  term  is  thus  2d  -2  and  the  last  term  has  the  order  Consequently  , 

we  require  2d  -2  - -<d,  or1*0  2/3. 

This  is  not  a necessary  condition  for  surface  stagnation,  for  if  we 
allow  x and  y to  be  of  different  orders,  we  may  have  no  cusp  at  all.  An 


example  of  this  latter  case  will  be  given  in  Chapter  V. 


i 


Chapter  IT 


In  this  chapter  we  develop  a representation  for  the  flow  under  a free 

^ / 

s+irface,  the  shape  of  which  is  a trochoid  (curtate  cycloid)-*  This  particular 
example  is  of  interest  because  it  is  related  to  an  early  attempt  to  find  the 
exact  equations  for  ocean  waves.  The  original  investigation  was  due  to 
Gerstner  [2]  and  was  later  duplicated  independently  by  Rankine  [?].  The 
main  difficulty  here  was  that  these  waves  were  rotational. 

If  one  attempts  to  solve  the  exact  equation  of  the  classical  ocean-wave 
problem,  keeping  only  terms  up  to  the  second  order,  then  one  is  led  to  the 
trochoidal  surface.  (Terms  of  the  first  order  yield  a sinusoidal  surface.) 
These  investigations  are  due  mainly  to  Stokes  [10]  and  Rayleigh  [6], 

The  approach  we  will  take  here  will  be  from  the  opposite  point  of  view 
with  respect  to  the  classical  theory.  Instead  of  specifying  that  the  bottom 
is  flat  and  then  searching  for  the  surface,  we  will  specify  that  the  surface 
is  a trochoid  and  find  out  what  kind  of  bottom  will  caus3  this.  It  is  known 
that  the  trochoid  does  not  represent  the  true  ocean  wave  form.  (This  chapter 
shows  this.)  So,  we  w-ni  answer  the  question,  "What  type  of  flow  does  produce 
the  trochoid  as  a free  surface  in  the  exact  irrotational  theory?" 

It  is  perhaps  helpful  to  have  a possible  physical  application  in  mind 
as  we  develop  the  theory.,  One  could  imagine  an  infinitely  long  washboard 
with  a sheet  of  water  flowing  along  it  at  a fixed  velocity.  The  theory  might 
apply  to  the  flow  of  water  in  a ^ver  if  it  were  sufficiently  straight  of 
~ / 

— ' ' Fritz  John  discovered  this  flow,  independently,  using  a method  similar  to 
that  of  this  paper.  His  general  method,  however,  treats  unsteady  flows  and 
this  example  was  a special  case  where  the  flow  was  made  steady.  John,  Fritz 
"Two  Dimensional  Potential  Flows  with  a Free  Boundary, « Communications  on 
Pure  and  Applied  Mathematics  t Vol.  VI  (1953)  pp°  1+97-503. 
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course,  uniform  m width  and  of  proper  speed,  depth,  and  bottom  shape. 

However,  the  main  practical  value  of  this  theory  is  its  possible  ap- 
plication to  the  study  of  ocean  waves  of  great  length  and  small  height, 
such  as  those  caused  by  underwater  seismic  disturbances.  As  is  show,  the 
roughness  of  the  bottom  becomes  infinitesimal  with  respect  to  the  roughness 
of  the  surface  if  the  ratio  of  wave  height  to  wave  length  is  sufficiently 
small.  The  flow  of  this  chapter  is  useful  more  as  an  indication  of  how  small 
this  ratio  must  be  in  order  to  have  the  existing  formulas  apply  than  as  a 
direct  source  of  expressions  useful  in  calculating  the  velocity  of  these 
waves,  for  velocity  expressions  have  existed  for  some  time  (Stokes)  which 
yield  good  results  for  waves  of  small  wave  height  to  wave  length  ratio. 

This  chapter  is  also  theoretically  valuable,  since  it  gives  an  estimation 
of  how  well  the  trcc'noidai  curve  approximates  the  profile  of  ocean  waves. 

The  approach  used  here  does  not  require  the  determination  of  a function 
A(<  ),  but  is  a direct  procedure  whereby  we  continue  the  flow  away  from  the 
given  free  surface.  Recall  that  in  the  existence  proof  we  assumed  that  a 
flow  was  given  together  with  its  free  s irface.  A consequence  was  that  the 
free  surface  was  analytic.  Now  it  is  knovm  that  an  analytic  function  is  com- 
pletely determined  if  it  is  known  along  a curve.  Thus  actually  th*  only  part 
of  a flow  that  needs  to  be  given  is  the  free  surface  itself,  for  then  by  means 
of  this  tneory  we  are  able  to  extend  the  flow  under  this  free  surface. 

There  are  several  metnods  applicable  to  this  problem  of  continuation; 
however,  they  are  all  variations  of  the  following  idea.  Since  the  free 
surface  is  analytic,  it  can  be  described  by  an  analytic  function  y = f(x). 

But  x = z - iy,  so  that  from  y * f(z  - iy)  we  can  solve  for  y to  obtain  y 


5' 
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% 2 

as  a function  of  z,  y - g(z).  If  we  multiply  equation  (?)  by  (~-)  wa  obtain 

1 * ai  TT"  (^>2 


whence 


— ^ - i ? a * 2±p.  53  . 

dz  i dz 

Replacing  A by  -g(z)  we  can  solve  this  differential  equation  for  - -^(z), 
which  is  the  required  flow  function.  Incidentally,  if  we  invert  — '■  -l>(z)  to 
obtain  z = z(~S)  t then  /(('£)  * - g(z(<)).  In  general,  to  get  from  the 
free  surface  to  the  actual  flow,  one  inversion  and  one  quadrature  are 
necessary. 

The  usual  representation  of  a trochoid  of  fixed  wave  length  2 If  , wave 
height  b,  and  of  variable  total  depression  c takes  the  form 

x = - Q ♦ b sin  6 , 

y * - c ♦ b cos  0 , 

where  we  require  c > b > 0.  The  condition  c > 0 follows  from  the  requirement 
that  the  free  surface  lies  below  the  x-sxis,  that  is,  A (~S)  > 0 for  "^real. 
The  definition  of  a trochoid  requires  c > b.  The  limiting  case  c * b = 1 
will  be  treated  in  the  next  chapter. 

Thus,  on  the  free  surface  we  must  have 

z = - 0 ♦ b sin  0 ♦ i (-c  ♦ b cos  0)  , 


where,  in  particular. 


= -(-c  ♦ b cos  0)  . 


We  note  first  that 
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"(f  ( ^ -c  * X)  d©  ■ (-1  ♦ b cos  ©)  dO  ■ -©  ♦ b sin  © ■ x. 


so  that  if  we  are  to  have 


on  the  free  surface,  we  must  have 


"J*  ^ ” -J'  (/^  -c  ♦ 1)  d©  . 


We  may  use  this  relation  to  detenoine  what  ^ “ .b  (9)  must  be  in  order  to 
yield  the  trochoidal  free  surface.  To  do  thi3,  we  first  differentiate  the 
equation  with  respect  to  ©,  obtaining 


/ 1 7?  „ 1 

1 • A d©  ^ 


c - 1 


or,  squaring. 


("57T  " ^,2)  (^")2  " ^ " 2(c’1)^  ♦ (c-D2 

Upon  replacing  by  & on  the  left-hand  side,  we  have 

- _) 

3^  (a^-)2  - <^)2 ' A2  -2(<=-i>  ♦ (=-D2 , 


or 


- 2 A j^2  - 2(c-l)^  ♦ (c-1)2  ♦ (%§-)2j  * 

We  now  insert  for  3 the  free -surface  function  c - b cos  0.  After  reduction 
(the  third  degree  terns  drop  out  by  virtue  of  the  Pythagorian  indentity)  we  have 


( 


^)2  - (2c-2b  cos  ©)  (b2  ♦ 1 -2b  cos  9)  . 


/ 


I! 

I 

1 
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In  general,  this  yields  a composite  elliptic  integral  for  £ . However,  if 

2 y 

we  choose  2c  - b ♦ 1,  the  right-hand  side  becomes  a perfect  square  and  o 

is  given  by  the  elementary  integral 


(b  ♦ 1 - 2b  cos  ©)  d©  . 


b?  ♦ 1 


This  special  choice  of  c corresponds  to  a choice  of  the  velocity  of  flow  along 

our  washboard.  It  is  quite  natural  to  expect  some  value  of  c to  be  favored, 

since  we  have  already  fixed  the  wave  length.  That  this  particular  value  of 

c is  the  ^best"  must  remain  a point  of  speculation.  The  numerical  work  of 

2 

tracing  stream  lines  in  the  case  of  2c  / b ♦ 1 becomes  unmanageable. 

Rough  calculations  gave  a strong  indication,  however,  that  the  value  c = - — ^ 
actually  did  yield  a smoother  bottom  than  other  values. 

In  order  to  have  Im(^5)  < 0 correspond  to  the  flow  region,  we  choose 
the  negative  sdgn  for  the  integral  and  obtain,  finally,  the  parametric 
representation  for  the  flow, 

5-  - (b2  ♦ 1)  © ♦ 2b  sin  © , 

b2  ♦ 1 

z •'-©♦bsin0*i  ( ?5 — ♦ t cos  G 

The  expression  for  the  velocity  is  particularly  simple  for  this  flow. 

We  have 


dj_ 

dz 


d2> 

~a©~ 

dz 

d© 


2 b cos  9 - (b  ♦ 1)  

-1  ♦ b cos  © - i sin  ?T 


,2  , / i©  . -i©N  . , 
b - b (e  ♦ e ) ♦ 1 


] - be 
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1 - be 


i© 
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whence,  since  the  velocity  is  the  complex  conjugate,  we  have 

_iQ 

velocity  = + i = 1 -be 

Since  a flow  must  be  represented  by  a single -valued  analytic  function 

* ~S(z)  t we  can  expect  trouble  in  the  above  flow  at  points  where  =?  is 
zero.  This  occurs  when 


, -ie  .. 

be  **  1 , 


or,  letting  6 


, when 


yj 

i ' (cos  ^ - i sin^")  " ^ . 


Since  b is  a positive  real  number,  we  require  £ * ?nTT  , and  the  trouble 
spots  are  hence  the  points 


6 - 2nTT  ♦ i 1q,t  _ . 

It  can  be  shown  that  the  flow  occurs  over  a multi-:,  belted  Hiepann  surface 
below  these  points.  This  is  indeed  unfortunate,  for  the  velocity  expression 
indicates  that  the  bottom  gets  smoother  as  we  go  deeper-.  However,  if  we 
are  to  remain  in  the  realm  of  physically  possible  flows,  we  must  stay  above 
the  images  of  these  points. 

The  actual  picture  of  the  flow  can  be  made  quite  easily.  7<"e  first. 

pick  a value  for  b,  where  2b  is  the  wave  height.  As  an  illustration  we 
1 


will  use  2b  * 


O' 


The  equations  become 


5 - I 0 ♦ -jkr  sin  Q , 


- © ♦ 


2 fT 


sin  0 ♦ i ♦ — y-TT  cos  0 ) . 


2 /f~2~ 


/ 


?! 

E 
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The  free  surface,  Im(^  ) » 0,  is  obtained  oy  taking  lm(0)  ■ 0 or  8 real. 
The  parametric  equations, 

1 


x - - © ♦ 


2{T 

y ’ 4 iff7" 


sin  0 


cos  0 


then  describe  the  free  surface.  A set  of  values  appears  in  Table  1.  The 

corresponding  curve  is  shown  in  Figure  1 on  page  31.  To  find  the  stream 

lines  below  this  free  surface,  we  seek  values  of  0 * £ 4 i^  whicn  will 

/ 

keep  the  imaginary  part  of  j a negative  constant.  Yfe  have 


<•  - l (E  4 i>3  ) ♦ 4= 

° ' ' *\ 2 


sin  ( £ 4 i > I ) 
' / 


“ - ^ * 4 sin«|  cosh  y - i(g  7j  - — — cos  £ sinh^/  ) J 

/ 1 2 


consequently  pairs  ( £ , 'Yj  ) must  be  found  such 


that 


9 1 

jr  7?  - — — cos  £ sinh?;  - k . 

V gT  5 / 


We  may  pick  for  k successive  values  only  up  to  k * .295,  fur  we  may  not  go 
beyond  0 *=  2nlT  4 i log  y , whence 

k < jr  log  2 fT — sinh  log  2 Y~2  ^ .295  • 

42 

In  Table  2,  we  compute  several  pairs  ( E ,yj)  for  k = .25  and  lie  resulting 
"bottom*  is  shown  in  Figure  1. 

In  this  particular  example  of  the  wave  motion,  the  "bottom"  is  rather 
shallow  and  quite  wavy,  hence,  not  very  useful  in  the  study  of  ocean  waves. 
However,  by  reducing  the  wave  height,  we  get  not  only  an  increase  in  depth  -■ 


i 


Figure  1 

The  trochoidal  free  surface  and  one  stream  line. 
.0^  units  abo/e  the  lowest  possible  stream  line. 


Bottom 

"wave" 

height 


Figure  2 

Bottom  roughness  and  maximum  depth  as  a function  of  v;ave  height. 
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that  is,  the  singularity  in  the  flow  moves  down  — but  percentagewise,  we 
get  a greater  smoothing  in  the  bottom.  To  compare  the  depth  of  singularity 
and  smoothness  of  bottom  for  various  values  of  wave  height,  b,  we  use  the 
relations  developed  above  v.ith  the  special  values  % -»  0 , '"j « - -log  b for 

the  crest  of  the  bottom  curve,  and  £ * TiT,  7j  **  ^ (where  the  value  of 

which  yields  the  same  imaginary  value  for  "t  as  * -lop  b)  for  the  trough 
of  the  bottom  curve.  The  quantity  i ♦ 1 - /J ^ is  then  the 

"wave  height"  of  the  bottom  curve.  Of  course,  2b  is  the  height  of  tne  surface 


wave  and 


<?/„  > - (iiii 


is  the  distance  from  trough  of  the  surface  wave  to  trough  of  bottom  curve, 
this  being  what  we  take  for  tne  depth.  It  is,  of  course,  not  the  usual 
definition  of  depth,  but  is,  we  think,  a useful  one  for  this  type  of  flow. 

A table  of  values  of  these  quantities,  together  v.ith  wave  height  to  wave 
length  ratio  and  percentage  smoothing,  is  found  in  Table  3 for  values  of 
wave  height  ranging  from  .fi  to  .0£,  the  wave  length  being  fixed  at  2 7/ . A 
composite  graph,  Figure  2,  displays  some  of  this  information.  Of  major  interest 
is  the  velocity  of  propagation  of  these  waves.  Since  in  the  treatment  the 
coordinate  system  moves  with  the  wave  crests,  we  must  look  for  the  mean 
velocity  of  the  water  past  the  observer.  Recall  that 


q * velocity  = ♦ i qj  “ ~ 


Thus,  the  absolute  value  of  the  velocity  along  the  free  surface  is 


(q  q)^^  * I i-2v  'os  9t  b4  (fer  0 real)  . 


The  velocity  at  the  crest  is  horizontal  and  of  magnitude  1 - b,  and  the 
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velocity  at  the  trough  is  also  horizontal  and  of  magnitude  1 •*■  b.  The  mean 
velocity  along  the  free  surface  is 


where  ^ is  the  tangential  velocity,  here  the  velocity  itself,  since  the 
.free  surface  is  a boundary^  and  c is  some  length  of  arc  rep-resenting  a 
complete  cycle.  Since  our  wave  is  symmetric,  we  may  use  a half-wave 
length.  Recall  that  on  the  free  surface 

x = 0 + b sin  0 , 

b*  * 1 . » 

y 2 * b cc£  0 > 

q = 1 - e'i0  » i - e-10  (for  0 real) 

= 1 - b cos  0 * ib  sin  0 * * iq.. 


Hence, 


f*!6*  * ^ dy  - - b cos  9)2  * b2  sin2  ©J  a© 

J y'dx2  ~ dy2  f{(l  - b cos  9) 2 * b2  sir.2  0 j 


-f 

Jo 


(1  r b - 2b  cos  e)  d0 


/"a . b‘  - 


ib  cos  9)  d0 


The  numerator  can  be  integrated  and  yields  the  value  - TT(  1 t b ).  Trie  de- 
nominator is  an  elliptic  integral  and,  through  the  substitution  0 = tr  - 2 y , 
takes  the  form 


l 


I 


2(l*b) 


V I 

:»)  \ > 1-  - 

lL  ( 


Ub 

(l*b)‘ 


sin 2<^  d<f  ■ 2(l+b)  E ( ^^— ) 


Thus  if  W2  take  the  negative  square  root  in  this  integral  so  that  the  signs 
will  cone  out.  right,  we  find  that  the  mean  velocity  of  the  fluid  on  the  free 


surface  is 


. -f  tB  <#>  r1  . 

m 2 (l*b)2  lb 


Recall  that  b is  the  amplitude  of  the  wave,  2b  is  the  wave  height.  It  should 
be  pointed  out  that  if  we  wish  to  interpret  this  flow  as  a progessive  wave, 
we  must  make  the  wavy  bottom  move  right  along  with  the  surface.  Consequently, 
■it.  is  better  to  think  of  this  flow  as  a river  flowing  over  an  undulating 
bottom.  If  the  bottom  of  a surface-wavs  flow  is  flat  or  infinite  in  depth, 
then  it  is  immaterial  whether  we  consider  the  waves  as  fixed  and  the  water 
in  motion,  or  the  water  fixed  and  the  waves  in  motion.  Thus  if  we  wish  to 
compare  the  flow  of  this  chapter  with  classical  ocean  waves,  it  is  best  to 
think  of  the  ocean  waves  as  progressing  up  a river  which  is  flowing  at  the 
velocity  of  propagation  of  the  waves.  The  bottom  of  the  river  must  bo 
and  the  river  must  be  quite  deep.  It  is  of  interest  to  compare  velocities 
obtained  from  the  various  existing  formulas  with  those  obtained  from  equation 
(13).  We  consider  these  formulas  for  fixed  wave  length  2^  and  for  g * 1- 
The  quantity  b is  the  amplitude  of  the  wave,  and  the  depth  is  assumed  to  be 
great  with  respect  to  the  wave  length  in  each  case.  From  the  linearized 

i.t.  . i.  l 1 4 X.  — A S'  ma  ne  4-  aa  i f*  f\  TT>r>X  T ^ S Tlfi  3.S  ©(3U3X 
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to  one,  Stokes'  first  approximation  (yielding  trochoidal  waves)  gave  the 
velocity  as  (l>b2)1//Z,  and  his  second  approximation,  carrying  terms  up  to  the 


i 


/ 


2 % h ’ 

fourth  order,  produced  the  relation  (1  t b ■*  b ) for  the  velocity. 
These  latter  expressions  are  compared  in  the  accompanying  table  with  values 
of  from  equation  (13)  for  various  wave  amplitudes. 
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(imb2) 

1.3013 

1.006 

1.020 

i.Olil 

1 

1.077 

1.12 

1.17 

j.  • 22  j 

J 1.25 

2 * h 1/2 

(lvb4  i-  4 b u) 

4 

1.0013 

1.005 

1.021 

1.01.9 

1.002 

1.15 

i.  • 2 j 

i.?ji 

l.li? 

Sn 

J..C02 

1 .00(3 

1.031 

] .066 

1.115 

1.18 

1.2-i  1 

L_J 

i . ? : 

1 .HO 

The  correspondence  in  velocities  is  quite  remarkable.  Fcr  1 rge  values 
cf  b,  however,  the  flow  of  this  chapter  satire t even  be  remotely  Compared 
with  deep-water  surface  waves,  sir-e.  +'>e  letter  is  very  wavy  ar.d  only  ; 
short  distance  bslcw  the  surface. 

The  results  cf  this  chapter  would  seem,  to  indicate  that  cither  there  is 
a remarkable  similarity  between  flows  over  a washboard  and  deep-ocsir.  waves, 
or  the  approximation s used  fcr  :os?.r,  waver  of  relatively  larye  -v ’plit.de  arc. 
no'  as  accurate  as  supposed. 
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Chapter  III 

Although  the  flow  described  in  this  chapter  car.  be  obtained  directly 
from  the  example  of  Chapter  II,  we  will  develop  it  in  all  detail  in  order  to 
illustrate  a different  method  of  attack.  When  one  inserts  a particular  function 
^ ( ~$)  into  the  expression  (8),  he  is  often  confronted  by  a highly  complicated 
integral.  In  order  to  force  this  integral  to  be  of  a simple  type,  one  could 
set  the  integrand  eoual  to  an  elementary  function.  In  this  chapter  we  will 
set  the  integrand  equal  to  unity. 

We  obtain,  then,  by  equation  (8),  the  system  of  equations 


(lU) 


1 , 


(15)  z ' - i/?  ♦J1  ♦ £ . 

Here  ^ serves  as  a parameter  to  connect  z and  . We  notice  that  on  the 
free  surface,  x is  equal  to  ^ . One  could  set  the  integrand  in  equation  (3) 

equal  to  any  function  of  /{ or  provided  the  conditions  stated  in  Chapter  I 

are  satisfied,  namely,  that  the  resulting  expression  obtained  by  solving  the 
differential  equation  for  ^ is  analytic  in  j,  positive  on  some  segment  of  the 

real  ^-axis,  and  satisfies  the  inequality  > 0 on  this  segment. 

If  we  solve  the  differential  equation  (lh)  for  we  obtain 


(16) 


A 


- l 


Equation  (15)  can  then  be  written 
(17) 


0 

r- 
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It  can  be  seen  from  equation  (16)  that  if  is  restricted  to  the  interval 
0 < ^ < 1/2,  is  real  and  varies  between  two  values  whose  difference  is  7TA* 
Thus  this  portion  of  the  real  -axis  will  correspond  to  the  free  surface.. 
Equation  (lU)  assures  us  that  > 0 on  the  frco-surface  segment. 


Analytic ity  is  obvious. 


The  integral  (16)  may  be  evaluated  by  means  of  the  substitution 
/{  - 1/2  cos^  9.  We  obtain 


i.. 


cos  t d9  = - 


9 sin  29 


and  from  (17), 


, 2 a 9 sin  29 

Z m - 1/2  cos  9 - ■£ J- 


where  these  equations  are  parametric  equations  for  the  flow.  With  this 
change  of  parameter,  the  free  surface  is  now  represented  by  real  values  of 
9,  and  the  interior  of  the  flow  is  obtained  by  allowing  o to  range  over 
certain  complex  domains  which  will  be  defined  later.  The  equations  may  be 
further  simplified  by  employing  the  identity  cos  Qm  1/2  (1  * cos  29)  and 
substituting^  for  29.  We  obtain 

f 

(18)  £ - - lA  ♦ siryf)  ; 

(19)  z - - 1/U  (/f*  sir yf)  - j-  (1  ♦ cos^) 

= - lA  1(1  * e’^)]  . 

The  shape  of  the  free  surface  may  be  found  by  seperating  z into  its  real 
and  imaginary  parts;  this  is  possible  on  the  free  surface  , sines  p is  real 
there.  We  obtain 


i 


/ 


1 
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(20)  x * - 1 A ( p + sin  £ < ) , 

y ■ - lA  ( l oos^3  ) . 

These  equations  are  easily  recognized  as  the  parametric  equations  for  a 
cycloid o 

We  must  now  carefully  determine  just  what  part  of  the  cycloid  represents 

the  free  surface.  Recall  that  /l is  restricted  to  the  range  0 < 1/2; 

consequently,  since  ^ ■ 1/2  cos^O,  this  corresponds  to  some  6 interval  of 

length  TT/2  which  we  may  take  to  be  the  interval  0 < 0 < JT/2,  or  in  terms  of 

tt.  For  this  range  of^5  , equations  (20)  will  describe  the 

half  cycle  of  the  cycloid  between  x «*  0 and  x ■ - 7T/U,  y varying  from  - 1/2 

to  0.  The  continuation  of  this  free  stream  line  from 0 to^  = - 7 T,  arises 

from  taking  the  negative  square  root  in  equation  (17).  The  segment  of  the 

real  A -axis  from  0 to  1/2  is  hence  traversed  twice.  This  corresponds  to 

the  segment  of  the  ^ -axis  lying  between  - TT/h  and  7T/u.  In  order  to 

investigate  the  flow  further,  we  look  for  the  image  of  the  remaining  part  of 

the  real -axis.  This  will  be  the  continuation,  as  a fixed  boundary,  of 

n 

the  free-surface  stream  line.  It  i3  no  longer  necessary  t o t ake /d  into 
consideration,  since  the  parameter p is  all  that  is  needed  to  connect  z and 
We  see  from  equation  (18)  that  the  continuation  ox  the  free  surface  as 
a fixed  boundary  corresponds  to  a curve  in  the  complex  -plane  along  which 
Im  (~£)  - Ira|-  lA  (^£>*  sin^  will  be  equal  to  zero  and  along  which 


Re 

•*  CD. 


j~  iMfi  * sin^)^  varies  between  - oo  and  - TT/h  and  between  TT/h  and 


Since  we  will  need  the  curves  for  which  Im  )■  constant,  we  will  develop 
these  curves  new, and  then  we  need  only  to  set  the  constant  equal  to  zero  in 
order  to  obtain  the  equations  for  the  continuation  of  the  free  surface.  If, 
in  equation  (18),  we  set^  • ^ ♦ i^7,  we 
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Im(^)  - Im  l/h  (£  * \7j  * si  cosh 7j+  i cos^sinh^) J 

- - l/li  ( yj  * cos  | sinh9^)  . 

Setting  our  constant  equal  to  - t:  w<s  obtain,  for  the  stream  lines,  the 
relation 


7j  ♦ cos  £ sinh  7jm  k 


whence 


(21)  cosf“tiHK?T  ' > ' ww“  vsinir 

/ 

In  order  to  find  the  image  of  these  lines  in  the  z-plane,  we  use  equation  (19) 


with  ft  m ^ ♦ i 7^  to  obtain 


. - . 77  -i  £ x 7 

X/4  15*  ly  ■*  mi  - e / e ' *J  * 


or 


— Uz  « £ * e^7  sin£  ♦ i(l  ♦?/  ♦ e^ccs  £ ) . 

To  write  x and  y as  functions  of  5 oirrfe  parameter,  we  may  use  equation  (21) 

1 jp 

together  with  the  identity  sin  £ * ♦ 1|l-cos  | , which  yields 

J“  * ' 

-hy  = 1 ♦ V * ^ * 

For  Yj  very  large,  we  see  that  x behaves  asymptotically  like  ety  and  that 

t 

y is  asymptotic  to  ly  . Consequently,  for  all  values  of  k.  including  k * 0, 
the  stream  lines  resemble,  asymptotically,  the  logarithmic  curve  y “ log  r- 


The  flow  is  illustrated  in  Figure  3.  Stagnation  points  occur  at  the  ends  of 
the  free  surface  where  the  cycloid  meets  the  fixed  boundary.  The  angle  here 


Tigure  3 

Th«  cycloidal  flow.  The  parts  of  the  line  k = 0 above  the  x-axis  represent  the  "roof" 
the  part  of  the  line  k = 0 below  the  x-axis  16  the  cycloidal  free  surface. 
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is  2TT/3.  There-  are  many  branches  tc  this  mapping,  but  the  only  one  involving 
a free  surface  and  locally  positive  pressure  is  the  one  shown  in  Figure  3* 

It  can  best  be  interpreted  as  r-  flr-v  f water  miner  i "roof"  shaped  like  the 
fixed  bourHarv . In  the  break  " be  rtr  jcf.w  the  water  drops  down  in  the  form 
of  a eye loid . 

We  find  the  velocity  by  ta>^ r.-.;  the  conjugate  of 


d z 

c/3 

~w 


= 1/2(1  - e1/5)  , 

1 v e-  V5 


w :ence 


Ir. 


1 t e'^c;^  - is"  ? sin  £ 


order  to  investigate  the  -.to.  hr  a.*,  .u  of  !.  ip.v  v,  out;  v>.d  : 
h power  ib-  s'.  p — • Wc  lavs 


-iiS  - 


p t sin ^5 


in  A - IT 


(/*-  ?T)3  (/S-  »T)  . 

— y ■■■  — • - “ 1 — ‘-y.  1 ■ + • • • . 


a ( XJ  w»\  2 / At  >i-\  3 

-U  id  T c-~V^)  . TT,  i Ti,  r 

I I */ 

;p*  A—  T7l  * ^ *.  ro  f/?a rrnti  — 

I/5  "l  ' 

-hr  - n - il  /3  ;:  .rr-)-  .•  :(£ ") 


and 


-i:;-  rr=  L^^ILl-  , vg')  . 


10' 


If  we  let  ^3  - fT‘  23*  j, 


then 


*/  By  f * 0(ot ) we  mean  that  hi  f k|oc!  , K “ cc  '-.start,  for  all  oC.  s iffisi 
near  tc  zero. 


• + i- 


3 i3«' 


~h<  _ rr=  r-fr—  + o(€ u)  . 

Now  the  free  surface  and  its  continuations  are  the  image  of  the  lines 
Im( $ ) = 0;  hence  fron 


0 - Ik(o)  * Tm(-h’fj-Tr)  = ^ sin  3©*  ♦ 0(£  l)  , 
we  see  that  necessarily  sin  39'  = 0,  or  9'  = 0,  * . 


Consequently,  we  have 


j«  - tt  = lidzpr 


0(6.  3) 


. r2  i29 1 
jTe  * 0l‘€  > 


_ r"  i(20»  + TT/2) 


t 0(6  ")  ; i - e~ 


Thus  in  the  neighborhood  of  ^9*  TT,  that  is,  when  £ — > 0,  the  argument  of 
-liz  -TT,  when  z is  tracing  out  a stream  line,  c..n  assume  the  values  TT/2.  7 TT/6, 
and  -Tt/6,  The  branch  corresponding  to  TT/2  is  the  cycloidal  free  surface,  and 
the  branch  corresponding  to  - TT/6  is  the  "roof."  The  third  branch,  cox-respond- 
ing to  ? TT/6,  yields  a flew  which,  for  position  pressures,  has  no  free  surface 
at  ail.  In  order  to  see  exactly  what  is  happening  rear  a branch  point,  it.  is 
sometimes  convenient  to  follow  a it  oar  by  stream  line. 

It  is  often  helpful  to  trace  out  various  stream  lines  in  the  parameter 
plane  by  composition  of  graphs.  Fbr  example,  in  tKis  flow  we  are  interested 
in  the  lines  for  which 


Ir.(-h  1'  ) - Y)  t sunh  7^c:s  j-  - constant  , ^ = ^ t i ~n  . 

Cne  plots  on  the  same  graph  the  11:  o:  7^  = constant  and  the  lines 

sir.h  Yj  cos  £ = constant  fer  equally  spaced  increments.  Ey  joining  the 

intersections  where  Yj  + si  nr.  Yj  cos  £ * constant,  we  construct  the  mages  < r 


i 


the  strean  lines.  Figure  h illustrates,  this  orocedure. 

For  a flow  to  be  physically  realisable,  vie  must  investigate  the  pressures 
involved  in  the  flow  to  make  certain  that  they  are  positive.  In  order  to 
do  this  we  consider  Bernoulli's  law  in  the  form 

constant  - ~ m l/2  q2  ♦ Im(z)  . 

If  the  pressure  is  to  remain  positive  as  vre  go  deeper  into  the  flow  region, 
the  quantity  l/2  q ♦ Im(z)  must  become  large  and  negative  there.  In  the 
case  of  this  example,  we  have 


1/2  q2  - 1/8  (1  ♦ 2e~7)  cos  £ ♦ e"2  ^ ) , 


Im(z)  » y ® -1./U  (1  ♦ ~SJ  * e ^ cos  £ ) . 

The  velocity  term  decreases  to  l/8  quite  rapidly  as  we  go  into  the  flow, 
since  vj  increases  there  (see  Figure  U)  . As  long  as  y i3  negative,  the  pressure 
is  positive,  but  if  y is  large  and  positive,  then  the  pressure  will  become 


negative.  The  latter  is  true  when  we  are  near  the  ••roof"  and  far  out  on 
either  side;  consequently,  this  flow  is  not  physically  realisable  in  the 
large , 
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Chapter  IV 

As  was  mentioned  in  the  introduction,  one  cannot  hope  to  pick  the  one 
function  X (^)  which  will  solve  the  ocean-wave  problem  unless  the  class  of 
functions  under  consideration  is  restricted  considerably.  In  this  chapter, 
we  will  consider  the  influence  that  the  problem  of  surface  waves  over  a deep 
ocean  has  on  the  selection  of  the  function  A ( i ) - First  we  will  show  that 
under  certain  assumptions,  ^ must  ba  a function  which  has  singularities. 
Then  the  exact  location  of  these  singularities  will  be  determined  and  the 
behavior  of  (if ) at  the  singularities  will  be  found  through  the  use  of 
conformal  mappings.  We  will  gain  enough  knowledge  about  ^ (^)  to  enable 
us  to  obtain  a good  approximation  to  the  solution  of  the  problem  of  waves 
over  an  infinitely  deep  ocean.  An  example  is  given  illustrating  the  results 
of  these  considerations. 

In  order  to  establish  our  theorems  about  the  behavior  of  we  must 

make  an  assumption  concerning  the  physical  flow.  We  assume  that  the  waves 
are  symmetric,  in  fact,  that  the  whole  flow  is  symmetric  about  the  perpendicular 
line  passed  through  the  crest  or  trough  of  the  wave.  When  one  considers 
the  fact  that  there  are  no  external  forces  acting  on  the  surface  — such  as 
wine --and  that  the  water  is  frictionless  for  the  purposes  of  this  theory, 
this  assumption  of  symmetry  seems  quite  valid.  A consequence  of  this  symmetry 
is  that  within  the  frame  work  of  this  theory , all  velocities  are  bounded 
and  bounded  away  from  zero.  This  does  not  mean  that  in  the  actual  ocean 
zero  velocities  are  not  possible,  for  it  must  be  kept  in  mind  that  our  frame 
of  reference  is  fixed  with  respect  to  the  waves,  not  with  respect  to  the  water 
at  great  depth.  Th’is  the  fact  that  the  velocity  of  the  water  is  bounded 
away  from  zero  means  that  in  an  actual  ocean,  the  velocity  of  the  water  never 
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attains  the  velocity  of  propagation  c £ the  wave.  This  assumption  appears, 
from  observation , to  be  valid  as  lore  as  the  rave  does  not  "break. ™ The 

breaking  of  a wave  apparently  occurs  when  the  velocity  of  the  water  at  the 

W 

crest  of  the  wave  exceeds  the  velocity  of  propagation  of  the  wave-*-' 

The  fact  that  the  velocity  is  bounded  has  the  sane  meaning  with  respect 
to  either  frane  of  reference. 

As  a result  of  our  assumption  of  sjmetry,  »e  can  show  not  only  that 
the  velocity  is  bounded  and  bounded  away  from  aero  in  the  entire  flow  region; 
we  can  determine  the  asymptotic  ievelopnent  of  z * 2.{"f  ) near  infinity  as 
well.  Suppose  we  consider  the  region  in  the  s-plane  bounded  above  by  a half 
cycle  of  the  free  surface  joining  the  trough  to  the  crest  of  the  wave,  and 
on  the  sides  by  the  two  sess-infinite  vertical  lines  dropped  from  the  trough 
ard  the  crest  of  the  wavs.  These  vertical  lines  are,  as  a consequence  of 
our  assumption  of  symmetry,  orthogonal  to  every  stream  line,  including  the 
free  surface. 

m ij 

Since  tie  flow  function  _5  * -i  (z)  is  an  anal.- tic  f=u.ctior.,  it  yields 
a conformal  sapping;  hence  the  region  defined  above  in  the  e- plane  naps  into 
a region  in  the  plane  wiu.cn  is  bounded  on  tep  by  the  image  of  the  free 
surface,  that  is,  by  a segment  of  the  real  I^-axis,  and  on  the  sides  by 
lines  that  are  everywhere  orthogonal  to  the  images  of  the  stream  lines.  But 
then  these  lines  imusfc  also  be  vertical  straight  lines  dropped  fron  the  ends 
of  the  segment  of  the  real  S -axis  that  represents  the  half  cycle  of  the 
free  surface  in  the  z-plane.  We  shall  call  the  region  in  the  s-plane,  and 
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— Stokes  has  estimated  that  a wave  will  break  when  the  ratio  of  wave  height 
to  wave  length  approaches  Ii7. 
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its  image  in  the  l^~plane,  a "ceil"  of  the  flow. 

We  will  locate  the  cell  in  the  z-plane  so  that  the  vertical  line  under 
the  trough  of  the  wave  coincides  with  the  negative  .imaginary  axis  from 
z *•  - ib,  b > 0,  toz  *=  - ioo,.  and  the  vertical  line  under  the  crest  of  the 
v/ave  runs  from  z “ air  - ia , 0 < a < b , toz  = d 7f-  i or>  „ In  the  plane , 
the  cell  will  be  bounded  by  the  segment,  of  the  real  "^-axis  joining  the 
origin  to  the  point  ^ a c7f  c > 0,  and  the  two  vertical  lines  dropped  from 
these  end  points.  The  origin  in  the  ^ -plane  shall  correspond  to  the  point 
z =•  - ib  in  the  z-plane.  Figure  5 on  page  62  illustrates  the  cell  as  it 
appears  in  the  z-  and  -J  -planes. 

We  now  map  the  cell  into  two  new  regions  by  the  mappings  w *=  e and 

v = e ■LJ/^o  In  the  v-plane,  the  cell  is  mapped  into  the  region  bounded  by 
the  semi-circle  of  radius  1 lying  above  tne  real  v-axis . and  the  segment  of 
the  real  v-axis  joining  y ■ - 1 to  v ■ ♦ 1.  The  image  of  the  cell  in  the 
w-plane  will  be  bounded  below  by  a segment  of  the  real  w-axis  joining 
w = a 1 to  i 5 b',  - 1 < a ' < 0 < b ' < 1 „ and  above  by  a smooth  arc  leaving 
the  real  w-axis  at  w = a'  vertically  upwards  and  arriving  at  w = b!  vertically 
downwards,  such  that  all  points  on  the  arc  are  less  than  unit  distance  from 
the  origin.  The  image  of  the  cell  in  the  v~  and  w-plan.es  is  illustrated  in 
Figure  6.  The  point  at  infinity  in  the  z-  and  ^ -planes  has  as  its  image  the 
origin  in  the  v-  and  w-planes.  The  cell  v/alls  correspond  to  the  segments 
of  the  real  axes  in  the  v-  and  w-planes. 

If  we  consider  the  domains  in  the  v-  and  w-planes  as  being  connected 
by  the  mapping  ->=  3>(z),  which  v;e  assume  as  known,  we  see  that  w as  a 
function  of  v is  real  on  the  real  axis  and  consequently  can  be  extended  into 
the  lower  half-plane  by  the  Schwarz  reflection  principle.  By  the  same  reasoning. 
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we  can  extend  v into  the  lower  half-plane.  These  extended  domains  are  the 
unit  circle  in  the  v-plane,  and  a region  bounded  by  a closed  smooth  arc  in 
the  w- plane  containing  the  origin,  but  contained  in  the  unit  circle.  Thus 
through  the  mapping  ^ * "^(z)  we  can  connect  these  extended  regions  in  the 
v-  and  w- planes  by  an  analytic  function  w(v)  which  is  achlicht  within  the 
image  of  the  cell,  and  transfbrms  the  origin  of  the  v-plane  into  the  origin 
of  the  w-plane,  and  the  positive  direction  of  the  real  axis  at  the  origin 
in  the  v-plane  into  the  positive  direction  of  the  real  axis  at  the  origin  in 
the  w-plane.  Since  the  mapping  is  schlicht,  the  derivative,  ^ , is  bounded 
and  bounded  away  from  zero  within  the  image  of  the  cell.  Furthermore,  the 
quantity  ^ is  also  bounded  and  bounded  away  from  zero  within  the  entire  image 
of  the  cell,  since  at  w = v = 0, 
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lim 

w-»0 


v dv  | 

W ' dwJwJsQ 


. (-)  al 

dv;  Jw=0  » 


dw 


and  we  know  that  ^ is  bounded  and  bounded  away  from  zero.  But  since 


1 dz 

— = 

o D 


d v dw 
c w dv 


it  is  clear  that  provided  o and  d are  finite  and  non-zero,  the  velocity,  q, 
of  the  flow  will  also  be  bounded  and  bounded  away  from  zero. 

The  mapping  function  w(v)  will  have  the  form 


w 


a«  v ♦ a?  v ^ , 


where  a..  is  real  and  positive,  and  the  series  is  convergent  for  1 v < 1.  If 
1 * ’ 


we  replace  w and  v by  e 


-iz/d 


and  e 


-i  S/i 


respectively,  we  obtain 


t 


-U9~ 


- iz/'d  _ -13/ c -125/  c -iS/c  M a2  -i5/c  N 

e - a, e + to..  a, e (if  — e 7 

x — x a^ 


Taking  the  logarithm  of  both  sides,  vre  find 

-1  I - log  (a.,  3-i  ^/c)  * log  (1  4 e'1*^0  4 ...)  , 


whence 


(22) 


z ” id  log  a. 


. d^  . ,da2  -i^/c 

t ♦ - o ♦ i— — e ■*  . 

L r*  a 


Finally,  we  have 


T 


dz  d d a2  ^ 
W7  * - r _ 5 e 


* . * - , 


which  shows  that  the  velocity  at  infinity  is  . 

•j3 — 

Since  velocity  - (^-) , we  may  conclude  that  in  the  equation 
(7)  z’  (z!  ♦ 2i  ^')  * , 

the  quantity  zl  « Vi?  is  also  bounded  and  bounded  away  from  zero  ir.  the 

entire  region  of  flow,  We  can  now  show  that  the  function  /?(^>)  must  have  a 
singularity  in  the  finite  part  of  the  flow  region  of  the  -^-plane.  We 
’•<111  first  show  thatj^j  < uj^J*  N in  this  region  for  suitable  constants 
M and  N.  If  we  solve  equation  (7)  for  we  obtain 


(23) 


d /?  1 


21 


From  the  form  of  this  equation,  and  the  bound  on  z”  from  above  and  below, 
we  see  that  the  right-hand  side  of  the  equation  is  bounded  by  some  constant 
M as  soon  as  A is  bounded  away  from  zero,  But  then 


/ 


Ql  „ C I i z’  L^l 

'vl  - j p^r  " 21  id-l 

< M j^j4,  N,  N = constant,, 

in  the  region  of  the  cell  excluding  3ome  circle  about  ^ » 0.  The  inequality 
is  obviously  satisfied  within  the  circle  about  /[  » 0.  Consequently > 
cannot  increase  any  faster  than  and  is  hence  finite  in  the  finite  part  of 
the  flow  region. 

Now  for  a flow  in  an  infinitely  deep  ocean,  the  entire  lower  half  of 
the  ^-plane  maps  into  the  flow  region,  and  the  entire  real  'S  -axis  represents 
the  free  surface.  Therefore  is  real  along  the  entire  real  I^-axis 

and  is  defined  and  bounded  by  M N in  the  lower  half  of  the  ^-plane. 

As  a consequence,  according  to  Schwarz's  reflection  principle,  we  define 
in  the  entire  upper  naif-plane  in  a unique  way  by  setting 
« "R  (^)  there.  Thus  ^ is  defined  and  bounded  by  M |^>j  + N in  the  entire 
^-plane.  If  we  were  to  assume  that  ^ (~$)  were  also  regular  in  the  entire 
^-plane,  then  by  an  extension  of  Liouville's  theorem,  ^ must  have  the  form 

A f ♦ B , A and  B constant. 


A T ■*  B 

j 


A and  B constant. 


But  this  choice  of  the  function  /(  (~S)  cannot  yield  a solution  to  the  surface- 

wave  problem,  for  it  generates  a free  surface  which  is  not  even  periodic. 

o 

Thus  we  must  conclude  that  /'V  cannot  be  regular  within  the  region  of  flow. 

We  will  discuss  the  type  and  location  of  the  singularities  of  /{(^> ) in  more 
detail  below. 

If  we  again  consider  equation  (?3)  and  use  the  fact  that  z'  is  bounded, 
we  see  that  if  0,  then  necessarily'  /{'  - oo,  that  is,  /l(~S)  is  singular 
at  the  image  of  this  point.  Thus  if  we  are  to  prevent  any  stagnation  points 


from  occurring  in  the  flow,  we  must  require  that 


p1- 


* 0 at 


/{  = 0, 
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\ Re  now  turn  our  attention  to  the  /(-plane  and  try  to  construct  the  image 

of  the  cell  there . The  part  of  the  boundary  that  corresponds  to  the  free 

surface  is  clearly  a segment  of  the  positive  real  axis,  for  we  require  /f  (if) 

to  be  real  and  positive  when  belongs  to  the  free-surface  segment.  To  find 

the  remain  ig  part  of  the  cell,  we  return  to  the  differential  equation  (7). 

As  we  move  along  one  of  the  vertical  boundaries  of  the  cell,  dz  and  d ~S  are 

clz 

both  pure  imaginary!  consequently,  along  one  of  these  lines  is  real. 

dPT 

This  implies,  in  effect,  that  the  velocity  is  real  along  the  vertical 

wall  of  a cell.  If  we  write  I for  known  imaginary  quantities  and  R-j,, 
for  known  real  quantities,  then  on  the  sides  of  the  cell  equation  (?) 
takes  the  form 


r (r  ♦ 2i  4jd_)  - Jl 


yjt  ’ 


2i 


or,  since  -j-  is  real,  we  have 


R3  ♦ Rk  d 


a - 1 


T7 


Thus  it  seems  that /r  is  real  along  the  walls  of  the  cell.  That  this  is  the 
case  is  a consequence  of  the  uniqueness  of  the  solution  of  equation  (7). 

If  one  starts  from  either  end  of  the  free-surface  segment  and  integrates 
down  the  side  of  a cell,  one  finds  that  in  determining  the  successive  approxima 
tions  of  Chapter  I,  page  lii,  all  quantities,  including  the  initial  values, 
are  real  so  that  necessarily  X is  also  real  along  these  vertical  walls  of  trie 


cell. 

We  can  show  that  our  cell  boundary  in  the  'A  — plane  is  partially  composed 
of  the  real  /?-axis  to  the  right  of  the  free-surface  segment,  and  the  real 
/^-axis  to  the  left  of  the  free-surface  segment  up  to  the  origin.  At  the 
origin  of  the  ^ -plane  we  must  have  a singularity,  for  as  we  found  in  the 
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If  we  solve  equation  (7)  for  A ' , we  obtain 


discussion  of  the  non-regularity  of  A , it  is  necessary  that  « 0 at  /?«  0. 


(2U) 


q.  - _5_L_ 

* Zi/T 


2z 


- /I)  . 


As  we  move  downwards  on  either  side  of  the  cell,  we  have  d-^  ” id  d 1~f/<  0, 
so  that,  since  z'  and  "Xare  positive  near  the  free-surface  segments,  the 


„ , Q 

sign  oi  u f ugicoo 


fith  that  of  ^ — — — — 


2z 1 


But  on  these  walls  of  the  cells,  the 


d ^ 1 2 

velocity,  q,is  real  and  equal  to  - — * i consequently,  the  sign  cf  d A agrees 

2 d z z 1 * 

with  that  of  /I  - . As  will  be  shown  immediately,  the  velocity  must  decrease 

as  we  move  down  the  cell  wall  under  the  trough  of  the  wave,  and  increase  as 

we  move  down  the  cell  wall  under  the  crest  of  the  wave.  At  the  trough  of  ths 

wave,  and  hence  d ^ are  zero  by  virtue  of  the  form  of  Bernoulli's  law  on 

12  Q 

the  free  surface,  namely,  -q  = - y,  and  by  the  fact  that  /I  ■ -y  on  the  free 

surface.  Thus  A is  stationary  there.  On  the  other  hand,  q decreases  as  we 

leave  the  trough  and  proceed  down  the  wall  of  the  cell.  Consequently, 

/\  — 1. _ becomes  positive  as  we  leave  the  trough.  Thus  d A is  positive  and 

A is  increasing.  But  if  q is  decreasing  and  /t  is  increasing  as  we  move  down 

the  wall  of  the  cell,  then  d A '"'ill  remain  positive  As  a result,  A must 

continue  to  increase  until  we  reach  a singularity.  The  right-hand  side  of 

equation  (2h)  is  regular  along  the  wall  of  the  cell  below  the  trough  of  the 

waveband  hence  there  can  be  no  singularity  in  A along  the  positive  real  axis 

to  the  right  of  the  free-surface  segment,  with  the  exception  of  the  point  at 

infinity.  Thus  the  wall  of  the  cell  below  the  trough  of  the  wave  maps  into 

the  real  axis  in  the  A -plane  lying  to  the  right  of  the  point  Am  b. 

Below  the  crest  of  the  wave,  just  the  reverse  holds.  Here  q is  increasing 

2 


so  that  the  quantity 


?-  -a 


will  change  from  zero  to  a negative  quantity  as 


i 
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vs  leave  the  crest  of  the  wave  and  move  downwards  on  the  wall  of  the  cell. 

Thus  d is  negative  and  /?  is  decreasing.  Again,  with  ^ decreasing  and  q 

increasing,  d^  remains  negative  and  consequently  /i  will  continue  to  decrease, 

at  a rate  bounded  away  from  zero,  until  we  encounter  a singularity  in  equation 

n 

(2k) „ This  singularity  occurs  when  f G.  Thus  the  segment  of  the  real  axis 
from  ^ * a to  * 0 corresponds  to  part  of  the  cell  wall  below  the  image  of 
the  crest  of  the  wave.  This  proves  that  /?<d)  has  a singularity  corresponding 
to  0 which  occurs  directly  below  the  image  of  the  crest  of  the  wave  in 
the  -plane.  Since  we  know  that  necessarily  “0  at  /(  0,  the  image 

of  the  cell  wall  in  the  /?-plane  must  turn  upwards  at  this  point.  We  are 
unable  to  find  the  exact  location  of  the  cell  wall  beyond  the  point  Am  0, 
however.  If  wo  could  find  it,  the  problem  of  ocean  waves  would  be  virtually 
solved. 

We  will  investigate  this  image  further  after  we  have  shown  that  q increases 
as  we  move  down  the  cell  wall  below  the  crest  of  the  wave,  and  decreases  as 
we  move  down  the  cell  wall  below  the  trough  of  the  wave.  To  show  this,  we 
must  assume  that  the  wave  rises  continuously  from  trough  to  crest.  This  is 
surely  the  case  for  a simple  wave.  If  we  consider  the  vertical  component  of 
velocity,  we  see  that  it  must  be  positive  on  the  free-surface  segment 

and,  as  a consequence  of  our  assumption  of  symmetry,  it  is  zero  along  the 
cell  walls  and  at  infinity.  But  P is  a harmonic  function,  and  consequently 
it  must  assume  its  minimum  on  the  boundary.  Thus  P > 0 within  the  cell. 
Furthermore,  if  a harmonic  function  assumes  its  minimum  within  a region,  then 
it  is  necessarily  constant.  Consequently,  (j>  is  actually  positive  within  the  cei3., 

If  we  move  infinitesimally  from  the  cell  wall  below  the  trough  of  the 
wave  into  the  interior  of  the  cell  region,  then' we  are  moving  in  the  direction 


i 

I- 


/ 


of  increasing  x and  increasing 


Consequently,  r > 0 on  this  cell  wall. 


Conversely,  if  on  the  cell  wall  below  the  crest  of  the  wave,  we  move  into 

the  interior  of  the  cell,  we  are  moving  in  the  direction  of  decreasing  x 

and  increasing  9 • Consequently,  ^ < 0 on  the  cell  wall  below  the  crest 

of  the  wave.  But  0 m 0 , so  that  0 > 0 on  the  cell  wall  below  the 

yx  1 xy’  l xy 

trough  and  < 0 on  the  cell  wall  below  the  crest.  This  is  just  our 

assertion  that  the  velocity,  which  is  equal  to  on  the  cell  walls,  is 

decreasing  as  we  move  down  the  cell  wall  below  the  trough,  and  increasing  as 

we  move  uunii  the  cell  wall  below  the  crest  of  the  wave,  since  when  we  are 

moving  downwards,  we  are  moving  in  the  direction  of  decreasing  y. 

Although  we  cannot  find  the  exact  image  of  the  cell  wall  beyond  A “ 0 

in  the  A -plane,  we  can  assert  that  it  must  lie  in  the  first  quadrant,  and 

must  be  an  analytic  arc.  Tie  will  show  that,  in  general,  /i  (~5 ) is  schlicht 

and  possesses  no  singularities  within  or  on  the  boundary  of  the  cell  besides 

'J*  -rT 

the  one  on  the  cell  wall  below  3 * c li  and  the  one  at  = co  which  have  already 
been  mentioned.  We  will  consider  'the  mapping  from  the  z-plane  to  the  /('-plane. 
Since  the  mapping  from  the  z-plane  to  the  3^-plane  is  regular  and  schlicht 
in  the  entire  region  of  flow,  what  we  assert  for  A ( z ) is  equally  true 
for  /)=  2(z(^))  = A(3), 

If  we  multiply  equation  (7)  by  the  quantity  we  obtain 

(25)  - 

d r 

On  the  right-hand  side  of  this  equation,  we  replace  by  - i (f>  . 

'n  — \ -) 

On  the  left-hand  side,  we  set  A*  A ^ 4 i A 2 ar‘d  dz  = dx  ♦ idy.  we  will 
consider  the  behavior  of  the  imaginary  part  of  equation  (2^)  as  we  move 
vertically  downwards  in  the  z— plane,  that  is,  for  dz  * idy,  dy  < 0.  Under 
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this  restriction,  the  imaginary  part  of  the  equation  is 


(26) 


^2  C'W--  % % ' 


ay 

We  have  already  seen  that  is  positive  in  the  interior  of  the  cell 

y 

and  on  the  free-surface  segment.  We  know  from  symmetry  that  Cp  is  zero  on 

y 

the  walls  of  the  cell;  consequently^  since  we  have  shown  that  the  velocity 

s 

is  bounded  and  bounded  away  from  zero,  must  have  the  3ame  sign  on  each 
wall  of  the  cell  as  it  has  at  infinity.  We  have  constructed  our  frame  of 
reference  so  that  ^ > 0 at  infinity,  and  hence  Px  is  positive  all  along 
each  wall  of  the  cell.  If  we  assume  that  the  free  surface  does  not  have  a 
vertical  tangent,  then  ^ is  also  positive  along  the  free-surface  segment. 
Consequently,  ^ is  positive  on  the  boundary  of  the  cell,  and  hence,  by  the 
minimum  principle  for  harmonic  functions,  x is  positive  throughout  the 
interior  of  the  cell.  Therefore  the  product  ? x ?y  is  positive  in  the 
interior  of  the  cell  and  on  the  free-surfacc  segment,  and  zero  along  the 
vertical  walls  of  the  cel] . 

Suppose  now  that  we  start  from  a point  on  the  free  surface  in  the  z-plane 
and  move  downwards.  We  wish  to  show  that  A-^  and  /<  p remain  positive  on  any 
such  vertical  line  dropped  from  the  free  surface  other  than  the  walls  of  the 
cell.  From  our  construction  of  the  image  of  the  C'sll  in  the  -plane,  we 
know  that  ? is  zero  on  the  free  surface  and  that  the  region  of  flow  lies 
above  the  free-surface  segment  in  the  /i-plane.  If  we  write  equation  (26) 
in  the  form 


(27) 


fy  ■ -k  ^2> 


we  see  that,  since  <ex  <py  is  positive  throughout  the  region  of  flow,  and 


/ 


-56- 


!|i 

?! 

i'j 


I' 

l-J 

't  ■ 


it 

I 

i 

? 


since  /1^  is  positive  on  the  free-surface  segment  and  /{ ^ becomes  positive 
as  we  leave  the  free-surface  segment  and  enter  the  interior  of  the  cell, 
the  quantity  ^ ^)  is  negative  as  vre  start  to  move  down  a vertical 

line  below  the  free  surface  in  the  z-plane.  Because  we  are  moving  in  the 
direction  of  decreasing  y,  this  means  that  the  product  ^2.^2  mus*  increasing* 
Now  > 0 on  the  free  surface,  /{  ^ - 0 on  the  free  surface;  consequently, 

^1^2  becomes  positive  as  we  leave  the  free  surface.  As  ?ong  as  A 2 is 
positive,  equation  (27)  shows  that  /{  g will  continue  to  increase.  Consequently, 
we  can  assert  that  ^ cannot  be  equal  to  zero  in  the  interior  of  the  cell, 
for  if  it  were,  then  the  quantity  “^2/^2  wou^  have  to  increase  to  zero 
through  positive  values,  an  impossibility.  With  /L,  positive, 
constantly  increasing  so  that  we  may  conclude  that  is  a1?n  positive 
throughout  the  interior  of  the  cell.  Thus  the  image  of  the  interior  of  the 
cell  lies  in  the  first  quadrant,  and  A is  non-zero  in  the  interior  of  the  cell. 

From  the  form  of  the  differential  equation(23) , we  see  that  /?can  have  a 
singularity  only  at  points  where  /(  * 0 or  at,  infinity.  We  have  seen  above 
that  cannot  be  zero  in  the  interior  of  the  cell'  consequently,  if  A is  also 
finite  there,  a fact  that  we  have  already  established  on  page  50,  then  /T 
must  be  regular  in  the  interior  of  the  cell. 

We  have  already  seen  that  /{  cannot  be  zero  on  the  wall  of  the  cell 

below  the  trough  of  the  wave.  We  also  found  that  A must  be  zero  at  a point 

on  the  wall  of  the  cell  below  the  crest  of  the  wave.  We  will  now  show  that 

/\  has  only  this  one  zero  on  this  cell  wall.  As  was  pointed  out  earlier, 

the  quantity  & is  zero  all  along  the  walls  of  the  soil.  Thus  if  we 

x y 

start  from  the  image  of  /?  «*  0 on  the  wall  below  the  crest  of  the  wave  and 
move  downwards  in  the  z-plane,  equation  (27)  assumes  the  form 


! 


l 
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ay 


^1  'V  " ■ 2 ^2  * 


We  have  already  established  the  fact  that  and  /<■  ^ are  positive  within 
the  cell.  Consequently,  on  the  boundary  of  the  cell  0 an<^  ^ 0* 

As  we  leave  the  image  of  /l”  0 in  the  z-plane  and  proceed  downward.  ^ 
must  assume  positive  values,  for  otherwise  the  cell  boundary  would  double 
back  on  itself,  and  since  the  interior  of  the  cell  must  lie  to  the  right  of 
the  boundary,  we  would  contradict  the  fact  that  A ± and  /?.  are  positive  in 
the  interior  of  the  cell.  But  if  become^  positive,  then  as  seen  before, 
the  quantity  must  be  increasing.  Consequently,  /i^  cannot  remain 

zero  as  we  move  away  from  the  image  of  /!=  0 in  the  z-plane,  but  must  also 
become  positive.  As  long  as  A p is  positive,  /?,/}  2 is  increasing.  Thus 
A cannot  be  zero  again  along  the  wall  of  the  cell  below  the  image  of  /?=  0. 

Even  more  is  true,  namely,  if  the  quantity  is  constantly  increasing, 

then  the  image  of  the  cell  wall  beyond  /\  = 0 must  intersect  every  hyperbola 
A * k from  left  to  right. 

Since  2 is  the  imaginary  part  of  , we  see  that  the  image  of  the 

2 

cell  wall  in  the  ^ -plane  is  always  rising.:.  Consequently,  in  the  mapping 
n o no 

from  z to  /i  or  from  1)  to  the  points  on  the  boundary  of  the  cell  map 

into  each  other  in  a one-to-one  continuous  manner.  But  if  this  is  the  case, 

#/  - 

the  mapping  is  schlicht  in  the  interior  of  the  cell  as  wellV  if  the  mappings 
from  z to  /?2  and  ^>to  AA  are  schlicht,  then  so  are  the  mappings  from  z to  A 
and  5 to  A . 

That  /?is  bounded  in  the  finite  part  of  the  ^-plane,  and  hence  can 


’v  i 


See,  fer  example,  Titcnmarsh,  The  Theory  of  Functions , Oxford,  1932,  p.  201. 
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have  no  poles  in  the  finite  part  of  the  ^ -plane,  is  an  immediate  consequence 
of  equation  (23),  as  was  pointed  out  earlier.  This  completes  the  proof  that 
the  function  2 c?)  can  have  only  one  singularity  in  the  finite  part  of  the 
cell  in  the  ^"-plane,  and  that  it  must  have  this  singularity  below  the  image 
of  the  crest  of  the  wave.  Further,  we  have  seen  that  arf)  is  a schlicht 
function  in  the  cell. 

We  have  already  shown  that  |;H4  N for  la~ge  values  of  /{ . We 

1 T? 

will  now  show  that  becomes  infinite  like  J at  infinity,  if  we  multiply 
equation  (2U)  by  2i/^,  we  obtain 


2i' /\  j\x  *»  i - z'  ( — 1 y - /{  ) . 


t 


3 . 


At  the  trough  of  the  wave,  we  have  seen  that  / i is  stationary  and  that  the 
velocity  decreases  at'  we  move  down  the  wall  of  the  cell  below  the  trough. 

As  we  move  down  this  cell  wall,  d"^  » id  Y dY<°,  so  that  is  increasing 


if y = -3_  < there.  As  ws  have  seen,  this  is  indeed  the  case.  If 

/)  2 is  increasing,  so  is  /{  , for  /{  is  real  and  positive  on  the  cell  wall 

n 2 

below  the  trough  of  the  wave.  With/i  increasing  and  --3-  decreasing  as  we 
move  down  the  cell  wall,  we  are  assured  of  finding  some  positive  constant 


such  that 


ii£ 

dS 


Z’  ( — =-jr 

2z’ 


ce  d^  ■ id  d < 0 , 


-/()<iiLfor/1>b*>b, 


> for 


b* 


Consequently, 


p\2  1 Mi|y|.  M2  , M2  = constant. 

Since  ""T^becomes  large  and  negative  as  we  move  down  the  cell  wall,  /? becomes 
large  and  positive  with  at  least  the  order  of  3 ' . If  we  write  equation 
(?)  in  the  form 


(28) 


(z  4 2i  ^ ) = r;r™  » 


wo  soo  t.hnt. 


Z ♦ 


2i/l  0 ( ^|"1/2)  d^  - 0 (|f|l/2)  . 

We  have  already  seen  by  equation  (22)  that  for  large  values  of  ^ , 

z » £ 0 (1)  ; 

consequently,  we  may  conclude  that 

2i^--|^*0  ( j^jl/2)  . 

By  reinserting  this  estimate  into  equation  (26),  we  have  that 


<z  ♦ 2i  Pk  ) - 0 (j^j  1)  , 


from  which  we  ma  ; conclude  that 


z ♦ 2i  ^ * O ( log  | i ) 


Thus  near  infinity 


/l-  | z ♦ 0 ( log 
i 

' ? 


j) 

i -[ 


1 (log  1 5 |) 


Thus  we  find  that  not  only  does  ^ become  infinite  at  infinity,  we  see  even 
the  bohavior  of  the  mapping  near  infinity.  Thewidth  of  the  cell  in  the 


/ 


* 
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/{ -plane  is  just  half  the  width  of  the  cell  in  the  z-plane  as  we  move  off 
to  infinity.  Since  the  width  of  the  cell  in  the  z-plane  is  one-half  wavelength, 
the  width  of  the  cell  at  infinity  in  the  A -plane  is  jj  wavelength. 

Finally,  we  shall  establish  the  behavior  of  Ml)  at  the  image  of  Am  0. 

We  find  that  we  must  have  a square  root  singularity  at  the  origin  in  the  A -plane, 
that  is,  the  image  of  the  cell  wall  turns  vertically  upwards  there.  We 
start  with  the  differential  equation  (23),  but  write  the  i ight-hand  side  as 
a single  fraction.  We  obtain 

d ?%  m 2iz'2  7 - i 
d " Uz ' ^ 

We  now  invert  this  differential  equation  and  consider  it  as  a differential 


- a.  j r» ^ .* sv 
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(29) 


Uz’  A 


A 


2iz'2/?  - i 


where  z'  is  considered  as  a known  function  of  the  dependent  variable  Jj>  . 
This  differential  equation  is  regular  at  = 0 and  consequently  will  have 
a regular  solution  there  of  the  form 


Q ♦ ^ * a2  ‘I2  ♦. . . 


— 0 at  ^ 0 


, whence 


From  equation  (29),  we  see  that  ^ 
has  a simple  zero,  so  a_  / 0.  Thus  near  /l  = 0 we  havei 


d ~5 

Also,  jj- 


* a 


which  proves  that  the  singularity  at  A m 0 is  a branch  point  of  the  first  order. 

This  concludes  our  investigation  of  the  function  Ml)  for  the  problem 
of  ocean  waves  over  an  infinitely  deep  ocean.  We  have  found  that  in  the 


. 


1 
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/(-plane,  the  cell  boundary  is  composed  of  the  real  ^-axis  and  a curve  which 
leaves  the  origin  of  the  /(-plane  vertically  upwards,  turns  immediately  to 
the  right  and  asymptotically  parallels  the  positive  real  /?- axis  at  a distance 
above  it  equal  to  one -half  the  width  of  the  cell  in  the  z-plane.  W -h  this 
information,  we  are  able  to  construct  a free-surface  flow  which,  because  oi 
our  lack  of  knowledge  of  the  image  of  the  cell  wall  beyond  /(»  0,  is  not  the 
solution  to  the  problem  of  ocean  waves,  but  which  does  yield  a good  approxi- 
mation to  the  solution  to  this  problem. 

We  are  concerned  only  with  the  mapping  connecting  the  ^ - and/? -plane, 
for  this  will  determine  our  choice  for  the  arbitrary  function  /?  ( £) . We 
assume  that  the  ends  of  the  free-surface  segment  in  the  ^ - and  /(-planes 
arc,  respectively >0  and  e 7T  , and  a and  b.  The  cell  boundary  ir.  these  two 
planes  is  illustrated  in  Figure  7,  page  62.  The  singularity  at  the  origin 
of  the  /^-plane  is  of  the  order  of 

In  order  to  determine  the  mapping  function,  it  is  more  convenient  to 
rotate  the  ^-plane  into  a w-plane,  90°  removed  from  the  plane,  by  the 
mapping  w * i ^ . Further,  it  is  more  convenient  to  place  the  singularity 
corresponding  to  A-  0 in  the  w-plane,  so  that  we  can  map  a polygon  in  the 
w-plane  into  the  real  axis  of  the  /(-plane.  Figure  8 then  illustrates  the 
figures  to  be  joined.  We  must  map  the  real  axis  of  the  /(-plane  into  a 
polygon  in  the  w-plane  thdt  is  composed  of  the  real  positive  w-axis,  the  part 
of  the  imaginary  axis  between  the  origin  and  the  point  ic  Y T,  the  line 
joining  icTT  to  a point  k ♦ icTT  , k > 0,  and  a ray  leaving  the  point  k ♦ icTT 
in  the  negative  real  direction.  We  will  require  the  vertices  of  t*"  polygon 
at  the  origin  and  icTT  to  correspond  to  two  positive  points,  b and  a,  b > a, 
on  the  real  /^-axis,  respectively » The  vertex  at  k ♦ ictr  will  correspo  id 


• SVK 


HftiT.  y --  . . r 


1 


-63- 


to  the  origin  in  the  A -plane.  A value  for  k cannot  be  assigned  — its  value 
is  determined  by  a and  t^and  c. 

Since  we  are  mapping  a polygonal  domain  onto  the  upper  half-plane,  wc 
may  employ  the  Schwarz-Christoffel  transformation.  The  angles  at  the  origin 
and  at  the  point  iciT  are  both  equal  to  YT/2.  The  angle  at  k ♦ icTT  is 
equal  to  -7 T.  Thus  the  Schwarz-Christoffel  transformation  takes  the  form 

* - if  /l  (P-a)~l/2  (/{■■  b)_l/2  d?  ♦ B, 

where  A and  B are  evaluated  by  the  point  correspondence  2-  b— -/w  « 0; 

^ * a — ->w  « icTT  . We  integrate  the  transformation  and  obtain 


Si)  (Pi-  b)  ♦ (a  ♦ b)  log  [ ) A - a*  ♦ \/\  - b*|  \ 

k ' J 

The  required  point  correspondence  yields 

A * ~-r  , B - - c log  (b  - a) 

u ’ ly 

whence,  returning  to  o b}  means  of  the  relation  - iw,  we  have 


♦ B 


(30) 


m - a'2»C'b  / (PI-  a)  ( A"  b)  - 2ic  log(  \Ja-  a + >/A-  b) 


ic  log(b-a) 


= - 1C 


1 a ♦ b 


b)'  » log  . 


This  last  equation,  or  rather  its  inverse,  is  our  choice  for  the  arbitrary 
function  /?(£). 

Wa  must  now  check  if  this  choice  of  the  function  P(^)  satisfies  the 
necessary  free-surface  conditions.  The  mapping  itself  takes  care  of  the 
requirement  that  A is  positive  for  ^ real.  That  2 on  the  free 

surface-segment  can  be  verified  directly.  We  have 


and,  consequently. 


(32)  y - <f^-)2  - - (a  - V*  sKZz.  b) 

1 2 

Thus  the  inequality  2X-  * becomes 

1 > (a  » b)2  (^-  a)  (3-  b) 

vr  - ^jt2 

or,  sines  0 < a < b on  the  free  surface,  we  have 

2 „ (a  ♦ b)2  (?-  a)  (b  -/?) 


The  right-hand  side  attains  its  maximum  for  A * 'fab"',  where  we  must  verify  that 

2 > (a  * b)2  ( /fab1  - a)(b  - -fat?)  c [(a  ♦ b)  ('Tb*  - j~a)]2 
2 -fab  2 

We  will  choose  a ♦ b » 1 and  c = 1 for  reasons  to  be  given  shortly.  Thus 

with  this  restriction  on  a and  b,  the  free-surface  condition  is  satisfied. 

Ar*  *3  2 

we  must  compute  ^ ) and  d^  as  functions 

of  /\  . The  first,  equation  (32),  has  already  been  found  in  connection  with 

d ^ a 

the  free-surface  inequality.  From  d ^ ~ d A and  equation  (31),  we  have 

d A 

d = 2clA  d A 

(a  ♦ b)  a)(  - b)’ 


Consequently, 


- lR  'Rr*  **  ^ 


- - i 


1 . P 1 . (a+b)2  (/^-a)(Pv-b)  ^ 2ci  "X  d/I  ^ 

a*b)  /O  -a)  (A  -b5 


^ X2  ♦ f 2c  A-  - (a*b)  ] /(  * ab 

(33)  (a^'JT7^n  d ^ • 

From  our  form  of  Bernoulli's  law  and  the  fact  that  /{  = - y on  the  free  surface, 
we  see  that  the  velocity  at  the  trough  of  the  wave  (i.e.,  for  /i  * b)  is  ^2b 
and  the  velocity  at  the  crest  of  the  wave  is  j2a.  If  we  compute  the 

quantity  ~q~  = and  take  its  limit  as  — >oo  along  the  real 


axis,  we  obtain 


— A 


a ♦ b 


Now  it  seems  quite  logical  that  if  our  flow  is  to  represent  a true  physical 
flow,  the  velocity  at  infinity  must  be  some  intermediate  value  between  the 
velocities  at  the  crest  and  the  trough  of  the  wave.  As  it  is  not  generally 
known  what  intermediate  value  should  be  taken,  we  will  let  the  wave  degenerate 
into  a flat  surface,  and  then  equate  the  flow  on  the  surface  with  the  flow 
at  infinity.  Thus  we  set  a = b and -obtain  the  relation 


If  we  set  c = 1,  this  is  satisfied  for  b = . Consequently, in  the  case 

c = 1,  if  there  is  no  wave  at  all,  the  velocities  at  infinity  and  at  the 
free  surface  are  equal  if  they  are  both  equal  to  one.  It  seems  logical, 
then, to  require  the  velocity  at  infinity  to  be  equal  to  one  when  there  are 
waves  on  the  surface.  Consequently, we  set  a ♦ b * I.  With  this  restriction, 


!• 
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ths  velocity  at  infinity  becomes  the  mean  square  of  the  velocities  at  the 
crest  and  trough  of  the  wave. 

We  are  now  ready  to  compile  our  formulas  and  compute  the  actual  flow. 

We  note  that  since  the  flow  region  is  above  the  real  A -axis,  the  arguments 
of  A-  b and  a change  from  zero  to  TT  as  we  pass  over  the  points  b and 
a,  respectively,  from  right,  to  left.  Consequently,  - b becomes  i Vb  - A 
as  we  descend  past  b,  and  n/a  - a becomes  i /a—  as  we  descend  past  a. 


Therefore , t.he 


the  mapping  in  the  z-piane  are 


and 


z 


P , I (x^a)(x+b) 

J V Tx-aJ  (Wj 


dx 


i A 


a 


for  0 < A < a . 


By  means  of  these  formulas,  we  may  compute  z for  values  of  A along  the  real 
A -axis.  This  corresponds  to  computing  z along  the  walls  of  the  cell,  up  to 
the  singularity  at  A = 0.  The  integrals  can  be  reduced  to  the  sum  of  an  ele- 
mentary integral,  an  incomplete  elliptic  integral  of  the  first  kind,  and  an 
incomp  ete  elliptic  integial  of  the  third  kind.  Unfortunately,  this  latter 
integral  cannot  be  found  in  tabulated  form.  It  was  necessary  to  compute  these 
integrals  by  machine  in  order  to  learn  anything  about  the  flow.  The  data  that 
were  computed  yielded  the  shape  of  the  free  surface  and  the  points  of  inter- 
section of  the  stream  lines  with  perpendiculars  dropped  from  the  wave  crest 
and  wave  trough. 
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To  find  the  points  of  intersection  of  the  stream  lines  with  the  vertical 
cell  walls  in  the  z-plane,  it  is  necessary  to  compute  pairs  of  values  of  X.  — 
one  less  than  a and  the  other  greater  than  b — which  yield  the  same  value 
for  the  imaginary  part  of  ^ . The  two  branches  of  the  .function  o = S'('A) 
(equation  (30))  to  be  used  in  this  calculation  are 


>/<i-  * log  ^ — > for  ^ > b. 


~ = Ir-  i{i4t  >/<* -*><*>  -a>  • i°e  — ~ - . 


for  0 < /(  < a. 


The  point  in  the  ...  -plane  that  corresponds  to  ^s0  is  hence 

As  the  wave  height  diminishes,  that  is.  as  b - a— > 0,  this  singularity  moves 

off  to  infinity  like  log  — - — . As  in  the  case  of  the  trochoidal  flow  of 

o — a 

Chapter  II,  the  flow  has,  in  general,  no  physical  reality  below  the  stream 
line  that  passes  through  the  image  of  'X  ~ 0 in  the  z-plane.  Since 

- 0 for  X * 0,  we  cannot  invert  z - z(“X)  into  a function  r\  - 'X(z) 
in  the  neighborhood  of  X ~ 0.  Thus  we  cannot  find  the  complex  potential 
i = - ( 'X  ( z ) ) in  this  region. 

It  is  interesting  to  note  that  the  inversion  of  :z(/^)  will  always  be 
impossible  at  X * 0 it'  we  are  to  require  that  all  velocities  are  bounded 
and  bounded  away  from  zero.  Recall  that  at  X ~ 0,  we  require  = 0 if 
the  velocity  is  to  be  bounded  away  from  zero.  In  order  to  be  able  to  in- 
vert z = z(X  ) at  a * 0,  we  must  have  -nr-  f 0 at  X = 0.  But  then 
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q 


dz 


o 


nhich  contradicts  the  fact  that  the  velocity  is  bounded  away  from  zero. 

In  the  example  of  this  chapter,  tha  velocity  at  <{  * 0 is  equal  tc  2 ifib. 

For  the  actual  computation,  a wave  height  of  .2  was  chosen,  so  that 
a * .U  and  b » .6.  The  results  of  the  numerical  integration  can  be  found 
in  Tables  5 and  6,  and  the  resulting  flow  is  illustrated  in  Figure  9 on 
page  69  e This  wave  is  a slightly  better  approximation  for  shallow  water 
waves  than  was  the  trochoidal  wave  of  Chapter  II.  In  the  case  of  the 
trochoidal  wave,  the  stream  line  through  the  singularity  at  A.  m 0 had 
a "wave  height"  of  .0585  units  when  the  surface  wave  had  a wave  height 
of  .2  units.  The  corresponding  figure  for  the  example  of  this  chapter 
is  .OUU  units.  The  wave  length  and  depth  are  approximately  the  sane  for 
the  two  flows. 

It  is  interesting  to  note  that  the  wave  heights  (A.  y of  Table  5) 
of  the  individual  stream  lines  do  not  decrease  monotonically  as  we  move 
down  away  from  the  free  surface,  but  exhibit  a small  oscillatory  behavior. 
This  might  indicate  that  the  velocity  of  the  wave  should  be  chosen  slightly 
differently.  On  the  other  hand,  it  might  ba  that  if  one  could  accentuate 
this  oscillatory  action,  one  could  obtain  a flat  bottom  at  the  minimum  of 
the  first  oscillation,  thus  obtaining  the  exact  form  of  waves  over  shallow 
water. 


If'  in  equation  (33)  we  were  to  choose  a value  of  c that  would  make 
the  numerator  a perfect  square  (c*  ■ ^ [ Va*  ♦ / b)  [a  ♦ b]),  then  the 
resulting  flow  is  the  trochoidal  flow  of  Chapter  II.  This  trochoidal  flow 
has  the  pro^rty  that  the  velocity  at  infinity  is  independent  of  the  wave 
height,  if  we  hold  the  wavelength  constant. 


7* 
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Figure  9 

The  cell  in  the  z-plane.  The  dashed  lines 
are  not  plotted  stream  lines  but  serve  only 
to  join  corresponding  known  intersections  of 
the  stream  lines  with  the  cell  walls. 


Chapter  V 


We  turn  cmr  attention  to  flows  that  are  non-periodic  in  nature  and 
which  have  vertical  fixed  boundaries  that  are  infinite  in  extent.  The 
vertical  fixed  boundaries  are  obtained  by  forcing  the  continuation  of 
the  free  surface  to  be  vertical.  If  we  again  consider  the  differential 
equation 

(7)  z (z  * + 2 iA  ’)  - 

we  can  determine,  to  some  extent,  what  type  of  function  A (S ) must  bs 
chosen.  Since  we  are  considering  a vertical  line,  dz  is  pure  imaginary, 
and  since  this  line  is  an  extension  o."  the  free-surface  boundary,  we  are 
on  the  real  if -axis,  and  d5  is  real.  Thus,  substituting  the  symbols 

R^j,  R^  for  real  quantities  and  I 2 for  pure  imaginary  quantities, 

equation  (7)  takes  the  form 

h (h  ♦ ^ > ■ rh  • 

or 

R2  ♦ R3dA  = 

on  the  vertical  fixed  boundary.  Consequently,  if  we  start  from  the  free 
surface,  where  A and  if  are  real,  and  integrate  along  the  real  vertical 
fixed  boundary,  the  solution  A (iS  ) must  be  real  by  virtue  of  the  fact 
that  the  coefficients  of  the  differential  equation  for  A are  real.  This 
is  the  only  possible  choice  of  A (^f ) , since  the  solution  of  equation 
(7)  with  prescribed  initial  data  is  unique.  Thus  we  must  choose'  real 
functions  for  A (if).  Furthermore,  if  we  consider  the  equation 

(3)  z - - iA  ^ -A  , 


we  see  that  the  vertical  fixed  boundary  extension  of  the  free  surface  will 


7 


1 ")  ' 

occur  when  -g-— — -a  is  negative,  and  will  be  infinite  in  extent 
i A j ip 

1 r\  / ^ *»  * ^ 

if  either/*  orj  y~  y *A  do  becomes  infinite  for  some  value 
of  5 on  the  extension  of  the  free  surface. 

From  the  standpoint  of  conformal  mapping,  we  see  that  we  must  require 
the  lower  half  of  the  S -plane  to  map  into  a portion  of  the  A -plane  that 
is  bounded  b-"  part  or  all  of  the  real  A -axis.  That  is,  the  real  o -axis 
must  map  into  part  or  all  of  the  real  A -axis.  If  there  are  nc  singular- 
ities on  the  free  surface  or  its  extension,  then  we  must  map  the  half-plane 
into  the  half- plane,  and  the  most  elementary  mapping  that  accomplishes  this 
is  the  mapping  A (S  ) " "S  . If  we  are  to  have'  a singularity  on  tho  real 
3 -axis,  it  should  correspond  to  A - 0,  since  the  differential  equation  (?) 
has  a singularity  at  this  point.  If  we  wish  the  entire  3 -axis  to  map  into 
a part  of  the  real  A -axis,  wa  must  have  a square  root  singularity  at  the 
origin  so  that  the  real  negative  3 -axis  is  doubled  back  onto  the  real 
positive  A -axis.  Thus  we  can  use  the  function  A (O  ) = ^ ^ • 

If  we  set  A m.~S  , equation  (8)  takes  the  simple  form 

z - _ i<  *\  J ~-~i  d$  . 


z * - i3  4J  h!  ~ 1 • 

A good  deal  of  information  about  the  flow  can  be  obtained  directly  from 


this  equation  without  going  into  the  computation.  We  see  that  for  large 
positive  _5  , the  integrand,  *j  — j,-  - 1,  is  pure  imaginary,  so  that  we 

are  on  a vertical  fixed  boundary.  As  ranges  from  1/ 2 to  0 we  niova 


along  the  free  surface.  Since  the  integral  involves  the  square  root  of 
the  inverse  of  3 , it  will  be  convergent  and  consequently  x will  move  a 
finite  aisuanoe  as  S traverses  the  interval  from  1/2  to  0.  During  this 
time,  y will  decrease  from  l/2  to  0.  When  3^  is  negative,  we  are  again 


. '»  . ?■  'Vit*  J-' 


/ 


i 


I 

i •*  i 


on  a fixed  boundary.  Thus  it  appears  that  the  flow  involves  two  separate 
vertical  fixed  boundaries  with  a free  surface  joining  them. 

As  the  computation  unfolds,  wo  obtain  first  a semi-fountain  in  two 
dimensions.  It  involves  a fixed  vertical  wall  with  a flow  which  moves 
up  the  wall,  departs  from  the  wall  at  right  angles  into  a cycloidal  free 
surface  which  meets  another  wall  tangentially,  and  finally  descends  down 
this  latter  wall.  This  flow  is  illustrated  in  Figure  10,  page  7U.  If 
this  flow  is  reflected  about  the  higher  wall,  we  obtain  a two-dimensional 
symmetric  fountain.  After  we  have  reflected,  we  may,  of  course, remove  the 
center  wall  so  that  the  free  surface  becomes  a complete  cycloid.  Figure  11 
illustrates  this  reflected  flow. 

It  was  mentioned  in  Chapter  I that  under  certain  circumstances  we 

could  have  a stagnation  point  on  a free  surface  where  the  characteristic 
2TT 

angle  of  — did  not  occur.  Figure  11  shows  an  example  of  such  a flow. 

The  center  of  the  cycloidal  free  surface  is  a stagnation  point  on  a free 
surface  where  the  free  surface  has  a continuously  turning  tangent,. 

By  choosing  another  branch  of  the  parameter  plane,  we  obtain  a flow 
which  occurs  in  the  region  of  the  z-plane  which  was  not  included  in  the 
above  flows.  This  exterior  flow  is,  in  fact,  the  analytic  continuation 
of  the  previous  interior  flow  across  the  free -surface  boundary.  It  comes 
down  uniformly  from  ♦ co  along  a vertical  fixed  boundary,  and  upon  meeting  j 

an  air  bubble  which  is  maintained  in  a slot,  separates  from  the  fixed 
boundary,  forming  a cycloidal  free  surface  with  the  air  bubble  and  finally 
flowing  off  to  - cd  along  a vertical  fixed  boundary  displaced  horizontally 
from  the  original  fixed  boundary.  ’ 


I 


/ 
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This  flow  is  illustrated  in  Figure  12.  Figure  13  is  the  reflection 
of  this  flow  with  the  center  wall  removed.  If  we  consider  ourselves  fixed 
with  respect  to  the  outer  flow  at  infinity,  we  obtain  a representation  of 
one  fluid  mass  rising  into  another.  The  free  surface  now  becomes  the 
interface  between  the  two  fluid  masses. 

This  latter  interpretation  might  be  cf  interest  to  meteorologists, 
for  one  finds  such  phenomena  in  the  atmosphere,  "when  an  air  mass,  which 
for  this  example  should  be  long  compared  to  its  width,  is  heated  near  the 
ground,  it  breaks  away  from  the  ground  and  rises  into  the  cool  surrounding 
air.  Of  0010*86,  in  the  actual  physical  occurrence,  a wake  is  formed  along 
the  outside  edges  of  the  rising  air.  However,  if  we  idealize  the  wake  by 
placing  partitions  in  the  air  which  rise  with  it  so  that  no  turbulence  is 
formed  on  the  sides,  then  the  occurrence  is  described  by  the  example  of 
this  chapter.  Figure  lU  illustrates  this  flow. 

Davis  and  Taylor  [l]  have  investigated  the  mechanics  of  large  bubbles 
rising  through  extended  liquids  using  an  approximate  theory.  Also  Taylor 
[12]  investigated  the  instability  of  liquid  surfaces  when  accelerated 
perpendicular  to  their  planes. 

To  develop  the  mathematics  of  the  flow,  we  return  to  the  equation 
developed  for  z, 

z - - i^  ♦ J dif , 

and  we  make  the  substitution  3 = -i  cos  9.  Therefore, 

■ - cos  9 sin  9 d9  , 
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so  that 


. 1 2 /* 

- l ^ cos  Q - 

i 2 . r . 

- j cos  S - I si 


tar.  9 cos  9 sin  9 d9 


sin  0 d9 


i 2 n Q sin  29 

- ~ 2 cos  9 - 2 * U * 


This  can  be  written  in  a neater  form  if  we  use  the  trigonometric  identity 
2cos^  0=1+  cos  20.  We  have  for  the  parametric  equations  of  the  flow. 

(3fc)  3 - iA(i  + cos  29) , 

z » lA(  - 29  + sin  29)  - j-(l  + cos  20). 

The  free  surface  is  obtained  by  taking  0 real  so  that  its  parametric 


equations  are 


X ■ l/h(-  29  + sin  20), 
y - - l/li(l  + cos  29). 

These  equations  are  easily  recognized  as  these  of  a cycloid. 

Going  into  the  complex  9-plane,  we  have,  setting  9 = £ + i ^ > 

^ = lA[l  + cos(2  \ + i 27p  ] 

■ l/li(l  + cos  2 £ ccsh  2rj  - i sin  2£  sinh  2 7)), 

so  that 

Im(  ) * - l/h  sin  2 ^ sinh  2 Yj  . 

The  stream  lines  are  thus  obtained  for 


or 


l/h  sin  2 sinh  21J  - c. 


Sin  2 § sinh  2 7)'  ? k = ^c- 


/ 
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On  the  free  surface  £ and  y must  satisfy  the  equation 

sin  2 j*  sinh  2 Tp  =0. 

The  actual  free  surface  is  obtained  by  taking  y = 0,  ij  arbitrary.  The 
continuation  of  the  free  surface  as  a fixed  boundary  is  obtained  by  taking 
2 £ * n VT  and  y arbitrary.  We  can  discover  where  the  images  of  these 
lines  are  upon  computing  z for  complex  values  of  0.  We  have 
z = - l/)j[i  ♦ i cos  29  + 29  - sin  29] 

= - l/ii [ 20  f i + i(cos  29  t i sin  29)] 

+ ssr\ 

- ~ l/h(2Q  tit  - e-^). 

Upon  setting  9 * £ t i , this  becomes 

z - - l/U[2  $ + 2±y  * i + i e 2l(*  + ) ] 

= - l/'uL2  £ + t i ♦ i e"  ' cos  2£  - e“  ' sin  2 £ j, 

so  that  we  have 

x = - l/h(2  $ - e-2^7  sin  2£), 
y = - l/h(l  + 2 y t e 2 7 cos  2£). 

We  see  that  the  fixed  boundary  continuations  of  the  free  surface  are  of  two 
types.  For  £ = n TT  and  y arbitrary,  the  fixed  boundaries  are  lines  drop- 
ping vertically  downwards  from  the  ends  of  the  cycloidal  domes.  For 
^ = n-^-  the  continuations  are  vertical  lines  passing  through  the  tops  of 
the  cycloidal  domes.  In  any  application,  only  one  dome  or  even  just  a 
half -dome  is  used;  however,  due  to  the  periodic  nature  of  the  equations, 
there  are  an  infinite  number  of  possible  domes  spaced  along  the  x-axis. 

Figure  15  on  page  77  shows  the  two  branches  of  the  parameter  plane  cf  0 
which  yield  the  internal  and  external  flows.  Table  7 contains  the  results  of 
the  computation  for  the  free  surface,  while  Table  8 contains  the  values  of  2 £ 
and  2y  together  with  the  corresponding  images  of  x and  y used  to  find  the 
stream  lines  k * i. 


***** 
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ve may  obtain  the  velocity  expression  from  the  relation 


From  equation  (3^)  we  have 

d "€  X 

= " ? sin  29' 

dz  1 1 i 

d©“  = “f  + 5COs29*2'  sin  2Gi 

whence 


cos  25  - 1 - i sin  29  2-2  cos  25 

2 2 — — 

If  we  use  the  identity  sin  25  * 1 - ccs  29  ■ (1  - cos  25) (1  -t  ccs  29),  this 

simplifies  to 

"q  * \ [sin  25  - i(l  tcos  26)]  = - ^ (1  t e1  2^). 

Finally,  inserting  for  25  the  quantity  2^  - i2  , we  have 

cf  = 2 e2  ^ sin  2 i;  - (I  + e2  ^ cos  2 | ). 

It  can  be  seen  from  this  expression  that  the  arrows  shown  on  the  stream 
lines  in  Figures  10  and  12  are  correct. 

As  was  mentioned  in  the  beginning  of  this  chapter,  another  possiole 
choice  for  the  arbitrary  function  /\(S)  is  the  quantity  With 

this  function,  we  obtain  a relation  which  represents  the  flow  past  a 
semi-infinite,  thin, vertical  obstacle.  On  the  down  stream  cr  right-hand 

| 

side  of  the  obstacle,  a free  surface  is  formed  which  leaves  the  obstacle 
tangentially,  but  quickly  turns  in  the  direction  of  flow.  On  the  upstream 
side  of  the  obstacle,  the  flow  parts  at  a stagnation  point  occuring  a short 


-sin  29 


sin  25  (1  - cos  29)  - j sin2  29 
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way  up  the  obstacle  with  the  upper  part  of  the  flow  pessing  off  to  infinity 
in  the  vertical  direction.  This  example,  which  is  shown  in  Figure  16  on 
page  80,  is  not  physically  significant  in  the  large,  for  the  free  surface 
continues  to  drop  as  we  go  away  from  the  obstacle,  so  that  very  large 
velocities  occur  at  points  far  removed  from  the  origin.  Furthermore , on 
the  upstream  side  of  the  obstruction  we  must  either  impose  a fixed  boundary, 
or  allow  negative  pressures  and  an  infinite  piling-up  of  the  water  in  front 
of  the  obstacle.  If,  however,  we  consider  only  the  local  properties  of  the 
mathematical  flow,  we  can  expect  to  obtain  a fairly  good  representation 
of  an  actual  physical  occurrence.  In  this  example,  the  best  physical 
picture  seems  to  be  the  high-speed  skimming  of  the  ocean  surface  by  a long, 
thin, vertical  obstacle.  Our  coordinate  system,  however,  has  been  chosen 
so  that,  mathematically,  we  are  riding  with  the  obstacle  and  the  ocean  is 
flowing  past  us. 

In  this  example  it  is  more  convenient  to  join  z and^  through  the 
parameter  A . From  X , we  have^  » A ^ , so  tha  ^ 

d5  - 2A  dA  , 

and 

d^_  = 1 
d^  - 2X“> 


whence 


( i-)2 


We  substitute  these  quantities  into  equation  (8)  and  obtain 

z - - iA  TT  ‘ ( T5T)2  2^ 

- - iA  ]/2A  -7  d A 

- - iA  ♦ (2A  - 1}3//2, 


><MW. 


Figure  17 

The  fourth  quadrant  yields  the  lover  part  of  the  flow. 


! 


The  behavior  of  the  free  stream  line  and  the  fixed  boundary 
extension  of  the  free  stream  line  are  now  apparent.  For  A.  > 1/2  we 
are  on  the  free  surface  which  extends  to  infinity  toward  the  right. 

For  "A  < 1/2  we  are  on  tne  fixed  vertical  boundary,  which  also  extends 
to  infinity,  but  does  so  vertically  upwards.  The  point  A = 0 is  a 
branch  point  of  the  mapping  and  we  can  expect  the  image  of  this  point 
on  the  vertical  fixed  boundary  to  be  a stagnation  point  where  the  stream 
lines  separate. 

The  exact  behavior  of  the  flow,  however,  can  be  determined  only 
through  calculation.  We  have  ^ ^ ^ > so  that  if  we  let  A.  ■ % ♦ i'p  , 

we  obtain,  for  the  images  of  the  stream  linos,  the  relation 

Im(j*  ) * 2^7/  = constant. 

These  are,  of  course,  equilateral  hyperbolas  with  the  coordinate  axes  as 
asymptotes,  as  illustrated  in  Figure  17.  If  we  set  the  constant  equal 
to  zero,  we  find  that  the  free  surface,  together  with  all  its  possible 
extensions,  is  the  image  of  the  ar.d  7^-axes. 

The  most  interesting  part  of  the  flow  is  obtained  by  considering 
the  fourth  quadrant  of  the  'A -plane.  As  has  already  been  pointed  out, 
the  actual  free  surface  is  the  image  of  that  part  of  the  real  A -axis 
for  which  ^ > l/2,  since  there  tne  free- surface  condition 

I 1^1 

-0’  or  ? 

is  satisfied.  The  free  surface  is  hence  given  parametrically  by  the 


equations 


, •>  3/2 


1/3  (2 { - 1)'/ 
> % - 1/2 


l 
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2/° 

Thus  the  free  surface  resembles  th«_  curve  y *=  x ' " . 

In  order  to  stay  within  the  region  of  flow  (the  fourth  quadrant), 
we  must  pass  below  the  branch  point  at  w 1/2;  consequently,  the  argument 
of  - 1 will  change  from  zero  t6  -77”.  Thus  the  argument  of  (2 A - l)"/r" 
changes  from  zero  to  - 3TT/2,  so  that 

(2  A - l)3/2  = i (1  - 2 A)3/2, 

for  0 < f\  < 1/2.  This  part  of  the  extension  of  the  free  surface  as  a 

fixed  boundary  is  therefore  the  segment  of  the  y-axis  from  y * - 1/2  to  y - l/3. 

The  point  C is  a branch  point  in  the  mapping.  If  we  were  to  con- 
tinue in  the  direction  of  the  negative  real  ^-axis,  we  v;ould  continue,  in 
the  z-plane,  to  rove  up  the  vertical  fixed  boundary.  However,  to  stay  in 
the  branch  of  the  mapping  within  the  fourth  quadrant  of  the  A -plane,  we 
must  turn  at  this  point  and  move  down  the  negative  77-axic.  We  set 
A = it2  and  obtain 

OX  O 7 /r' 

2 » - t ♦ ^ (1  ♦ 2it“~;  *■ . 

This  is  the  stream  line  that  meets  the  vertical  fixed  boundary  at  a 
stagnation  point.  A3  this  particular  line  is  of  interest,  we  will  obtain 

its  parametric  equations.  Let  us  cali  it  the  “zex'o,e  stream  line. 

0 io< 

The  argument  of  1 * 2 it"  = r e is  between  zero  and  JT/2  and  is 
2 

given  by  cA  = arctan  2t  . In  terms  of  cC , r and  t,  we  nay  write 


- - t2  .4  ( 


3 • • 

cos  ~ 00  1-  1 sin  -|C<)  , 


whence 


t2  - ] r3/2  sin  3 * 


1 3/2 

-r  r ' 


cos  ^ °<  . 


/ 


Employing  the  identities 


cos  1 oc=  (2CC3«._  1)  (— 4°-g-— )l/?,  (0<oc<  TO, 
sin  | « = (2cos  ♦ 1)  (—  (o<ac<lT) 


and  the  relations 


we  obtain 


cos  <x. 


1 ,2  1/2  t v 1/2  / ^ . v 

= » t - 7 (r  - 1)  ' , (r  > 1 ) 


i-i 


- A <r2  - l)1/2  - A (2  . 1)  (i iI)V2. 


» - § r3/2  <|  - 1)  (--^>1/2,  (r  » 1) 


or,  upon  simplification. 


x *=  - (?-^)l/2  [(i-p)1/2  ♦ ^ (2  ♦ r)  ], 

y = ^ (2  - r)  i1  £ r)"'/2,  (r  > 1). 


The  derivative  assumes  the  simole  form 


dy  „ v r - 1 

dx  vT  ♦ ^T-rr 


This  further  extension  of  the  free  surface  accordingly  leaves  the  y-axis 
at  right  angles  in  the  negative  x direction,  turns  immediately  downwards 
and  finally  ha3  asymptotically  the  slope 

If  we  look  for  the  flow  above  thi3  "zero"  stream  line  we  must  consider 
the  third  quadrant  in  the  *X -plane.  The  images  c>f  the  stream  lines  here 
correspond  to  choosing  positive  constants  for  the  imaginary  part  of  : 
consequently,  we  will  soor.  run  into  negative  pressures  if  we  stray  far  from 


the  "zero"  stream  line 


/ 


« • I 

* I 


it 
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Cons  idering  again  the  lines  Im(^)  = 0,  but  dwelling  now  in  the 

third  quadrant  of  the  A -plane,  we  first  retrace  our  steps  up  the  7p-ax is, 

In  the  z-plane  we  are  of  course  coming  back  up  th  "zero"  stream  line, 

only  now  we  are  on  top  of  it.  TThen  we  reach  t)  origin  in  the  /\  -plane, 

2 

w»  turn  and  run  out  the  negative  real  axis.  Thus  7}*  Q,t  ” - P and 

we  have  0 9 -> /o 

3 " ip“  ♦ T i (i  + 2p  j ' . 

rm.  J x i x . r»  x ' r «.  - t , _ __  n /h 

io  oiic;  pal  b Ox  ox;c  y — d uuvc  y E x/  _>  • 

The  picture  of  the  flew  now  becomes  clear.  The  fluid  flows  up  from 
the  left  at  an  angle  of  approximately  h5°,  strikes  a plate  obstacle  which 
separates  the  flow  into  a part  which  goes  upward,  and  a part  which  flows 
downward  and  soon  separates  from  the  plate  into  a free-surface  wake.  The 
stagnation  point  where  the  flow  divides  is  the  point  where  the  "zero* 
stream  line  meets  the  y-axis. 

Turning  to  the  velocity,  we  use 

dj  **  2A  d/i  , 

dz  * (~  i ♦ '/z?  - 1)  dA 

to  obtain 


a - a ♦ ia  „ - 
H H1  M e.  Ndz  / 


2T 


i ♦ -J2  - 1 

..  i * »£>-  i1 


On  the  free  s’irface, 


q-;  - 


O **  — "1  « 


J 


% > V2 
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On  the  vertical  segment  from  y ■ - 1/2  to  y * 1/3, 

" 0> 

q2  “ - 1 * /l  - 2?  , 0 < J < 1/2. 

Finally  along  the  '•zero"  streai.  line  we  have 

n - fr  - 1>V2 

Qt  ' O ' > 


- - 1 ♦ (?L4-^)1/2, 


where  this  r is  equal  to  the  r used  in  the  parametric  representation 
of  the  "zero"  strcar:  line.  Along  the  vertical  segment  from  y ■ l/3  to 
y ” ♦ oo  we  can  use  the  relations  developed  for  the  lower  segment. 


1 U»«  1 WMW  ir  w «,«.n  .4 
oauv>|/u  uiia  u i iv  *•  1'uuj.u  a. 


~~  ~ — A 4-~ 

x*  '-'in  v-/ 


Table  1 


/ 


Points  on  the  free  surface  for  the  trochoidal  flow. 


y - T5  * 


cos  £ 


£ 

-X 

-y 

0 

0 

.210 

.10 

.065 

.213 

.20 

.129 

.21? 

.50 

.330 

.252 

.75 

.510 

.305 

1.00 

.703 

.372 

1.25 

.9lU 

.U5o 

rr/2 

1.22 

.563 

1.75 

l.Uo 

.627 

2.00 

1.66 

.711 

2.25 

1.97 

.786 

2.50 

2.29 

.8U6 

2.75 

2.62 

.888 

2.90 

2.62 

.906 

3.00 

2.95 

.913 

7T 

3.lU 

.917 

3.30 

3.36 

.91.3 

3. 50 

3.62 

.895 

3.75 

3.95 

.853 

U.oo 

U.27 

.793 

t.5o 

U.85 

.637 

5.oo 

5.3t 

• U60 

5.5o 

5.75 

.312 

2TT 

6.28 

.210 

Table  2 


Pciints  on  the  'stream  line  k 

* .25  for  the  trochoidal 

flow. 

-1  (l.l25?-  .25).  . fc 

♦ — e y sin>  ; y m - V - 

'i.rr 

c.  re 

• rr  ♦ 

1 

003  1 .V67  sinhw  r X 5 

• j 

T5 

7 

I 

-X 

-y 

.1363 

rr 

rr 

1.10U 

.1365 

3.05 

3.01 U 

1.103 

.137 

2.98 

2.92 

1.100 

.lU 

2.76 

2.61 

1.080 

.15 

2.U1; 

2.18 

1.026 

.17 

2.06 

1.71 

.936 

.2 

1.7U5 

1.32 

.839 

.3 

1.15 

.7lU 

.669 

.81 

,U2 

.599 

..5 

.56 

.25 

.568 

,6 

.335 

.123 

.55U 

.65 

.215 

.06U 

.551 

.67 

.16 

* 

0 

VA 

O 

.550 

.69 

.07 

.021 

r-’  1 . n 

- 

.695 

0 

0 

.5i*9 

e"?  cos£ 


Table  3 


Values  of  7^  and  7^  together  with  Ay,  2b,  £ smoothing  and 
wave  height;  wave  length  ratio. 


■ - log  *r.? 

~ ♦gb_ 

v?  72  _ 
b * i 

sinhl?  0 * 

f *- 

V- , * 5 *y 

•' 1 b2  ♦ 1 

- be 

V2 

» 

4l'i  ♦ 1 - 

■h 

smoothing 

2b  -Ay 

2B 

100;  wave 

length  = 2TT. 

b 

°?1 

% 

A y 

2b 

% 

smoothing 

ratio 

.U 

3 

C\ 

• 

.1136 

.6U57 

.8 

19.3 

1:8 

.3 

1.20U 

.237U 

.Ul38 

.6 

31.0 

1:10 

.2 

1.609 

.U90 

.2075 

.U 

U8.1 

1:16 

.15 

1.897 

.7lU 

.1233 

.3 

59.0 

1:21 

.1 

2.303 

1.07 

.0585 

.2 

70 

1:31 

.05 

2-996 

1.73 

.0160 

.1 

8U 

1:62 

.025 

3.6889 

2.lil35 

.00393 

.05 

92 

1:125 

i 


Table  a 


K>i 

i 


Depth 

2 * 

~ yb  “ *•  y value  of  bottom  below  trough  ■ 4 — “ 4 

. b2  ♦ 1 

- y “ - y value  of  surface  at  trough  ” — * + b 

£J  ^ 

depth  - ys  - yb 


b 

"yb 

-ya 

depth 

r 

.U 

1.1U17 

.98 

.1617 

s 

.3 

1.1628 

3 

CO 

* 

.3179 

%£ 

JBfc*  ' 

.2 

1.336 

.72 

.616 

- 

.15 

1.532 

.661 

.871 

.605 

«r?  * • 

a*'*?  B 

.1 

1.867 

1.262 

V«_' 

k 

ft 

.05 

2.513 

.551 

1.962 

■ 

.025 

3.193 

.525 

2.663 

i! 


‘-fcMtffcv  ****** 
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Table  5 


Intersections  of  the  stream  lines  with  the  walls  of  the  cell. 


f .h 

A > -6 

-y  trough 

-y  crest 

L y 

,h 

.6 

C 

.6 

.U 

• c. 

.395 

.60219 

.2508 

.831 

.677 

.15U 

.38 

.6080 

.U897 

I.0U9 

.936 

.113 

.36 

.6158 

.6711 

1.223 

1.128 

.095 

,3U 

.622U 

.7972 

1.3k3 

1.260 

.083 

,32 

.6280 

.8936 

1.U31 

1.360 

.071 

,30 

.6335 

.9706 

1.515 

1.UU0 

.075 

28 

.6385 

1.0335 

1.593 

1.505 

.088 

26 

.6U20 

1.086U 

I.62U 

1.563 

.061 

,2k 

.6U50 

1.129U 

1.653 

1.60U 

.OU? 

20 

.6510 

1.1971 

1.728 

1.673 

.055 

16 

.6560 

1.2iiii6 

1.783 

1.722 

.061 

12 

.656U 

1.2772 

1.798 

1.755 

.0U3 

08 

.6602 

1.2979 

1.827 

1 7TA 

.051 

oli 

.6607 

1.3C92 

1.832 

1.788 

.(Mk 

0 

.6611 

1.3126 

1.835 

1.791 

.olUt 
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Table  6 


feirs  of  values  on  the  free  surface., 


■yr 

-y 

0 

.6 

.3U8 

.5  9$ 

.70  T 

.58 

l.ooU 

.56 

1.2U7 

.5U 

1.U62 

.52 

1.66^ 

.50 

1.863 

,U8 

2.066 

M 

2.285 

.uu 

2.5'U7 

.1j2 

2.8U2 

.li05 

3.127 

.h 

Tahla  7 


Jbints  on  a cycloid. 

Ux  - - 20  * sin  20 
Uy  “ - 1 - cos  20 
0 < 20  <TT 


20 

UiX 

-hj 

§5 

1 0 

0 

2.0 

! -3 

.ooU 

1.96 

•5 

.02 

1.8b 

1.0 

.16 

1.5U 

£ 

1.57 

.57 

1.0 

Rw 

2.lU 

1.30 

.U6 

~ 

1 

2.6U 

2.16 

.12 

5 

2..8U 

2.$h 

.oU 

■fr' 

r 

3.lii 

3.1U 

0 

Tablo  8 


?.£  and  2 ^ 3uch  that  sin  2£  sinh  2^  ■ ^ and  corresponding 

Ux  and  Ly. 

l*x  « - 2f  ♦ e~2y  sir  2$  ; Uy  « - 1 - 2J?  - e~2?  cos  2| 


2$ 

2 2? 

i*x 

l*y 

3.1U 

♦ 00 

-3.il* 

- 00 

3.oU 

2.31 

-3.03 

-3.21 

2.8U 

1.20 

-2.76 

-2.01 

2.6U 

.91 

-2.1*5 

-1.56 

2.1U 

.56 

-1.66 

-1.25 

1.67 

.1*9 

-1.06 

-1.1*3 

1.57 

.1*8 

- .95 

-1.1*8 

1.U7 

.1*9 

- .36 

-1.55 

1.0 

.56 

- .52 

-1.87 

.5 

.91 

H 

• 

1 

-2.26 

.3 

1.28 

- .22 

-2.55 

.1 

2.31 

- .09 

-3.1*1 

-.X 

-2.31 

- .91 

-8.7l* 

-.3 

-1.28 

- .78 

-3.18 

-.5 

-.91 

- .69 

-2.27 

- 1.0 

-.56 

- .1*7 

-1.39 

-1.1*7 

-.1*9 

- .15 

-.67 

-1.57 

CO 

• 

i 

- .05 

-.52 

-1.67 

-.1*9 

.05 

-.35 

~2.ll* 

-.56 

.61 

C"? 

-2.61* 

-.91 

1.1*5 

2.09 

-2.81* 

-1.28 

1.76 

3.71* 

-3.0U 

-2.31 

2.03 

11.36 
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